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Abstract. The semiclassical mechanics of the Wigner 6j-symbol is examined 
from the standpoint of WKB theory for multidimensional, integrable systems, 
to explore the geometrical issues surrounding the Ponzano-Regge formula. The 
relations among the methods of Roberts and others for deriving the Ponzano- 
Regge formula are discussed, and a new approach, based on the recoupling of 
four angular momenta, is presented. A generalization of the Yutsis-type of spin 
network is developed for this purpose. Special attention is devoted to symplectic 
reduction, the reduced phase space of the 6j-symbol (the 2-sphere of Kapovich 
and Millson), and the reduction of Poisson bracket expressions for semiclassical 
amplitudes. General principles for the semiclassical study of arbitrary spin 
networks are laid down; some of these were used in our recent derivation of the 
asymptotic formula for the Wigner 9j-symbol. 



PACS numbers: 03.65.Sq, 02.20.Qs, 02.30.1k, 02.40.Yy 

1. Introduction 

The Wigner 6j-symbol (or Racah VF-coefiicient) is a central object in angular 
momentum theory, with many applications in atomic, molecular and nuclear physics. 
These usually involve the recoupling of three angular momenta, that is, the 6j- 
symbol contains the unitary matrix elements of the transformation connecting the 
two bases that arise when three angular momenta are added in two different ways. 
Such applications and the definition of the 6j-symbol based on them are described by 
Edmonds (1960). More recently the 6j- and other 3nj-symbols have found applications 
in quantum computing (Marzuoli and Rasetti 2005) and in algorithms for molecular 
scattering calculations (De Fazio et al 2003, Anderson and Aquilanti 2006), which 
make use of their connection with discrete orthogonal polynomials (Aquilanti et al 
1995, 2001a,b). 

The 6j-symbol is an example of a spin network, a graphical representation for 
contractions between tensors that occur in angular momentum theory. The graphical 
notation has been developed by Yutsis et al (1962), El Baz and Castel (1972), Lindgren 
and Morrison (1986), Varshalovich et al (1981), Stedman (1990), Danos and Fano 
(1998), Wormer and Paldus (2006), Balcar and Lovesey (2009) and others. The 6j- 
symbol is the simplest, nontrivial, closed spin network (one that represents a rotational 
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invariant). Spin networks are important in lattice QCD and in loop quantum gravity 
where they provide a gauge-invariant basis for the field. Applications in quantum 
gravity are described by Rovelli and Smolin (1995), Baez (1996), Carlip (1998), Barrett 
and Crane (1998), Regge and Williams (2000), Rovelli (2004) and Thiemann (2007), 
among others. 

Alongside the Yutsis school of graphical notation and the Clebsch-Gordan school 
of algebraic manipulation there is a third approach to the evaluation of rotational 
(SU(2)) invariants. The third method, sometimes called chromatic evaluation, grew 
out of Penrose's doctoral work on the graphical representation of tensors and is closely 
related to knot theory. We will not have further occasion to mention this school, see 
Penrose (1971) for its introduction, Rovelli (2004) for an overview and Kauffman and 
Lins (1994) for its full development. 

The asymptotics of spin networks and especially the 6j-symbol has played an 
important role in many areas. By "asymptotics" we refer to the asymptotic expansion 
for the spin network when all j's are large, equivalent to a semiclassical approximation 
since large j is equivalent to small ft. The asymptotic expression for the 6j-symbol 
(the leading term in the asymptotic series) was first obtained by Ponzano and Regge 
(1968), or, more precisely, they obtained several formulas, valid inside and outside 
the classically allowed region and in the neighborhood of the caustics. In the same 
paper those authors gave the first spin foam model (a discretized path integral) 
for quantum gravity. The formula of Ponzano and Regge is notable for its high 
symmetry and the manner in which it is related to the geometry of a tetrahedron 
in three-dimensional space. It is also remarkable because the phase of the asymptotic 
expression is identical to the Einstein-Hilbert action for three-dimensional gravity 
integrated over a tetrahedron, in Regge's (1961) simplicial approximation to general 
relativity. The semiclassical limit of the 6j-symbol thus plays a crucial role in simplicial 
approaches to the quantization of the gravitational field. 

For all these reasons, the asymptotic formula of Ponzano and Regge for the 
6j-symbol has attracted a great deal of attention. Ponzano and Regge obtained 
their formula by inspired guesswork, supporting their conclusion both with numerical 
evidence and arguments of consistency and plausibility. The formula itself is of the 
one-dimensional WKB-typc, a reflection of the fact that the 6j-symbol fundamentally 
represents a dynamical system of one degree of freedom. 

The Ponzano- Regge formula was first derived by Neville (1971), using the 
recursion relations satisfied by the 6j-symbol and a discrete version of WKB theory. 
Similar techniques were later used by Schulten and Gordon (1975a, b), who also 
presented stable algorithms for evaluating the 6j-symbol numerically. A proof of 
a different sort was later given by Biedenharn and Louck (1981), based on showing 
that the Ponzano- Regge formula satisfies a set of defining properties of the 6j-symbol. 

More recently there have appeared more geometrical treatments of the 
asymptotics of the 6j-symbol, that is, those based on geometric quantization (Kirillov 
1976, Guillemin and Sternberg 1977, Woodhouse 1991), symplectic geometry and 
symplectic and Poisson reduction (Abraham and Marsden 1978, Arnold 1989, Marsden 
and Ratiu 1999) and other techniques. Among these are the works by Roberts (1999) 
and by Charles (2008). In addition, the 6j-symbol has been taken as a test case 
for asymptotic studies of amplitudes that occur in quantum gravity (Barrett and 
Steele 2003, Freidel and Louapre 2003), in which the authors developed integral 
representations for the 6j-symbol as integrals over products of the group manifold. 
There have also been quite a few other studies of asymptotics of particular spin 
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networks, including Barrett and Williams (1999), Bacz et al (2002), Rovelli and 
Speziale (2006), Hackett and Speziale (2007), Conrady and Freidel (2008), Alesci 
et al (2008), Barrett et al (2009), among others. We also mention the works of Gurau 
(2008), which applies standard asymptotic techniques (Stirling's approximation, etc) 
directly to Racah's sum for the 6j-symbol; of Ragni et al (2010) on the computation 
of 6j-symbols and on the asymptotics of the 6j-symbol when some quantum numbers 
are large and others small; and of Littlejohn and Yu (2009) on uniform approximations 
for the 6j-symbol. 

In addition there has been some work on the q-deformed 6j-symbol, important 
for the regularization of the Ponzano-Regge spin-foam model (Turaev and Viro 1992, 
Ooguri 1992a, b) and for its possible connection to quantum gravity with cosmological 
constant. In particular, Taylor and Woodward (2004) applied the recursion and 
WKB method of Schulten and Gordon to the g-deformed 6j-symbol. The results 
are geometrically interesting (the tetrahedron of Ponzano and Regge is moved from 
M 3 to S 3 when the g-deformation is turned on), but it seems that at present there is 
no geometrical treatment of the asymptotics of the g-deformed 6j-symbol, analogous 
to what is available for the ordinary 6j-symbol. There is also the recent work of Van 
der Veen (2010) on the asymptotics of general g-deformed spin networks, which treats 
the problem from the standpoint of knot theory and representation theory. Among 
other things, this work creates a broad generalization of the Schwinger-Bargmann 
generating function of the 6j-symbol. 

In Aquilanti et al (2007) we applied multidimensional WKB theory for integrable 
systems to the asymptotics of the 3j-symbol, and in this paper we apply similar 
techniques to the 6j-symbol. These methods bear the closest relationship to the 
works of Roberts (1999) and of Charles (2008). The point of this paper is not another 
derivation of the Ponzano-Regge formula, although one is provided, but rather to 
clarify the relationship among some of the methods used in the past, to reveal useful 
calculational techniques, and to lay the basis for the development of new results. 
Among the latter we mention our own work on uniform approximations for the 6j- 
symbol (Littlejohn and Yu, 2009), our recent derivation of the asymptotic form for 
the 9j-symbol (Haggard and Littlejohn, 2010), both of which relied on techniques 
explained in this paper, and our work on the Bohr-Sommerfeld quantization of the 
volume operator in loop quantum gravity (Bianchi and Haggard, 2011). Previous and 
current work on the volume operator includes Chakrabarti (1964), Levy-Leblond and 
Levy-Nahas (1965), Lewandowski (1996), Major and Seifert (2001), Carbone et al 
(2002), Neville (2006), Brunnemann and Rideout (2008, 2010) and Ding and Rovelli 
(2010). 

In addition, this paper is distinguished by its use of what we call the "4j-model" 
for the 6j-symbol, in contrast to the "12j-model" used by Roberts (1999). The 4j- 
model is less symmetrical than the 12-model, but it is closer to the manner in which 
the 6j-symbol is commonly used in recoupling theory. In addition, in loop quantum 
gravity (Rovelli 2004) angular momenta represent area vectors, which in the case of 
four-valcnt nodes correspond by Minkowski's (1897) theorem to a tetrahedron. In 
this context the 4j-model is closer to the applications than the 12j-model, indeed, it 
played an important role in the work of Bianchi and Haggard (2011). 

In this paper we refer to Aquilanti et al (2007) as I, for example writing eqn. (1.13) 
for an equation from that paper. We note two errata in I, namely, a(x) in (1.89) should 
read a(u), and jz and j'4 should be swapped in the 6j-symbol in (1.112). 
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2. Spin network notation 

We begin by explaining our notation for spin networks, which is based on that of Yutsis 
et al (1962) with modifications due to Stedman (1990). At the end of this section we 
compare our conventions for spin networks with others in the Yutsis tradition. 

2.1. The 3j -symbol and Wigner intertwiner 



mi 



Jl J2 J3 
mi mi ms 




J3 m-3 



Figure 1. The 3j'-symbol contains the components of the standard three- valent 
intertwiner. 

The 3j-symbol is a number that can be regarded as the components of an 
intertwiner W : Cj 1 (& Cj 2 (g) Cj 3 —> C with respect to the standard basis \jimi) ® 
\j2m2) 55 IJ3TO3), as indicated in Fig. [TJ In this paper Cj denotes a carrier space for 
unitary irrep j of SU(2), so that dimCj = 2j+l. We call W the "Wigner" intertwiner. 
We will gradually explain the features of Fig. QJ as we proceed. 

The standard notation for the 3j-symbol is on the left of Fig. [1] while the central 
diagram is the standard Yutsis spin network for the 3j-symbol, with small arrows 
presented as in the Yutsis notation. The indices (mi, r»2, m;j) in the central diagram 
are covariant, that is, they transform under rotations as the components of a dual 
vector (in contrast to an ordinary vector). In a Hilbert space we regard ordinary 
wave functions or ket vectors as "vectors," while bra vectors are regarded as "dual 
vectors." Thus, contravariant indices are those that transform as the components of 
a vector. In the Yutsis notation the arrows indicate the transformation properties of 
the corresponding m index, and there are rules for "raising and lowering" indices, that 
is, reversing the direction of the arrow. The rules do not, however, make use of the 
metric as in ordinary tensor analysis. Our definition of the arrow (explained below) 
is different from that of Yutsis, but designed so that the two notations agree as much 
as possible. In particular, our notation for the 3j-symbol is the same as the Yutsis 
diagram in Fig. [TJ 

A trivalent node of a spin network such as those illustrated in Fig. [TJ is assumed 
to have a positive or counterclockwise orientation, unless otherwise indicated (thus we 
dispense with the + sign used by Yutsis). 



2.2. Bras, kets and scalar products 

The diagram on the right of Fig. [TJ makes use of the standard basis vectors \ jimi) in 
Cj i , i = 1,2,3. The spin networks for these basis vectors and their duals are shown in 
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\jm) = ^ m (j m \ = ~) m* 

Figure 2. Spin networks for bras and kets. 

Fig. [5] The large, broadly open arrow is a "chevron" (Stcdman 1990). When pointing 
outward (inward), the chevron indicates a ket (bra) vector. Also shown in Fig.[3]is the 
spin network for an arbitrary vector in Cj, and the dual bra vector (tp\ obtained 
by Hermitian conjugation. In the Dirac notation it is customary to label bras by 
the same symbol as kets, it being understood that the two are related by Hermitian 
conjugation. This convention is so deeply ingrained that we dare not change it. But in 
spin networks there are two different ways of converting kets into bras and vice versa, 
and this presents some notational challenges. One can see in Fig. [5] that Hermitian 
conjugation applied to bras and kets is notationally the changing of bra chevrons to 
ket chevrons and vice versa, and the starring of identifying symbols, with a double 
star being removed. Full rules for Hermitian conjugation of any spin network are given 
in Sec. 

The lines of a spin network will be referred to as "edges," including cases like 
those shown in Fig. O 

An edge of a spin network ending in an unstarred m index represents a contraction 
with the basis ket \jm), so the m index transforms under rotations as a covariant 
index. The explicit insertion of basis kets may be seen in Fig. [TJ An edge ending in 
a starred m index represents the insertion of a basis bra (jm\, so starred indices are 
contravariant. Thus, the star on an m index indicates its transformation property, and 
invariant contractions can only take place between a pair of starred and unstarred m 
indices. 

m = y^—r = r^—< 

Figure 3. Orientation of a spin network does not matter. 

As illustrated in Fig. [3l the orientation of a spin network on the page does not 
affect its value. The spin network in the figure has been rotated by 180°. 

(</#> = r 3 « 3 $ = 0*—^ v 
(</#) = (r -2- — v)* = r 

Figure 4. Spin networks for scalar products or contractions. Illustration of rules 
for complex conjugation. 

As illustrated in the final diagram in Fig. [T] and in Fig. 21 when a bra chevron 
and a ket chevron are juxtaposed, it represents the scalar product or contraction. In 
effect, the bra chevron acts as a receptacle for the ket chevron, and vice versa. After 
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the contraction the two edges may be joined, with a small arrow remaining to indicate 
which was the bra and which the ket in the contraction. One might suppose that the 
star would carry the same information, but, as shown below, it is possible to change 
the direction of the arrow without changing the stars. Figure HI also presents another 
example of Hermitian conjugation (complex conjugation, in this case). Since the small 
arrow represents the contraction of a bra and a ket chevron, when the chevrons are 
reversed, the direction of the arrow changes. 

{jm\jm') = m* — <^ — - — m' = m* — i — m' 

j 

= S mm > = ml* — m 



Figure 5. Orthonormality relations for basis vectors. 



As a special case of the scalar product, the orthonormality relations of the basis 
vectors are illustrated in Fig. [5] The final spin network follows from the symmetry of 
the Kronecker delta, or, alternatively, by the reality of 5 mm i and the rules for complex 
conjugation. 

The j labels on the edges of the spin networks indicate which carrier space Cj 
the ket or bra lies in, or in which carrier space a contraction has taken place. If 
there are distinct carrier spaces with the same j label, then additional distinguishing 
information must be supplied. 



2.3. Intertwiners 




Figure 6. Spin network for the Wigner or 3j'-intertwiner. 

An SU(2) intertwiner (the only kind we are interested in) is a linear map between 
two vector spaces that commutes with the action of SU(2) on the two spaces. The 
only case of interest here is where the target space is C, consisting of scalars, that is, 
invariants under rotations. The Wigner intertwiner W : Cj 1 ® C J2 <Ei Cj 3 —> C is of this 
type. But a linear map from a Hilbert space to C can be thought of as a dual or bra 
vector, for example, we can associate the map W with a bra vector (W\ belonging to 
C* ± ®C* 2 <8>C* 3 • The spin network notation for (W\ is shown in Fig. [SI The components 
of the Wigner intertwiner, that is, the 3j-symbol, are obtained by inserting basis kets, 
the ket chevron first, into the bra chevrons of the intertwiner, as in the final diagram 
of Fig. [U 

More generally, a C-valued intertwiner on a Hilbert space H can be regarded as 
an 5i7(2)-invariant bra vector on this space, that is, a member of H*. By Hermitian 
conjugation we obtain an SU (2)-invariant ket vector in %. Thus, there is a one-to-one 
correspondence between the subspace of H. of rotationally invariant vectors and the 
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set of C- valued intertwiners on W. The subspace Z introduced in Sec. 13.41 below is a 
subspace of this type, consisting of rotationally invariant vectors. 

2.4- Tensor products and resolution of identity 

J2\jm}{jm\ = < 3 rnm*^L < = ^ 1 ( 

m m 

Figure 7. Resolution of identity, and an illustration of the outer product. 

The outer product of a ket with a bra is represented in spin network language 
simply by placing the spin networks for the ket and the bra on the same page, as 
illustrated in Fig. [7] The orientations of the bra and ket spin networks is immaterial, 
but in the figure they have been placed with their m indices adjacent in order 
to emphasize the summation (a contraction over two indices, one covariant, one 
contravariant). The final diagram is the spin network for the identity operator on 
space Cj . 



7 7 7 
.4 — mm* — -L = -L. 



Figure 8. Summing on m to join two edges. 

Figure [7] illustrates another technique, the replacement of a sum on m by a joining 
of edges. The general usage is shown in Fig. [8] The directions of the arrows and the 
stars on one of the m's must be coordinated as shown for the identity to be used as 
shown. The small arrow in the final diagram in Fig. [3] (a fragment of a spin network) is 
a reminder of the directions of the arrows before the sum. This small arrow is omitted 
in the final diagram of Fig. [7J (the identity diagram) because the chevrons already 
indicate the direction of the edge. Recall that edges are also joined on contracting a 
bra with a ket, as in Fig. [4] 

< 3 « J * = <^—^ 
^*^_ ^ 1 </ = 0* — I ^ 

r^-^< — <f^V> = 4>*^— i> 



Figure 9. Summing on m to join two edges. 

The use of the identity diagram is illustrated in Fig. [9] The first row illustrates 
its action on kets (as a map : Cj Cj), where bra and ket chevons are combined as in 
Fig- HI The small arrow is omitted after the joining of the bra and ket chevrons because 
the remaining chevron indicates the direction of the edge. The second row of the figure 
illustrates its action on bras (the map : C* —> C*), and the third row illustrates its 
action as a map : C*®Cj — > C, that is, the scalar product map : (4>\ <g> \ip) >-> {</)\ip}. All 
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of these usages are encompassed by the same spin network. The identity spin network 
can also be seen as an element of C,- <E> C*. that is, the space dual to C* <g> C,- on which 
it, viewed as a C-valued linear operator, acts (the third line of Fig. [9]). 

In general a spin network has some number of edges that terminate in incoming 
or outgoing chevrons. Examples are the identity diagram in Fig. [7] and the Wigner 
intertwiner in Fig.|6l In all cases there are multiple interpretations of the spin network 
as a linear operator mapping one vector space to another, depending on how many 
of the incoming and outgoing chevrons have ket and bra chevrons plugged into them 
(specifying the domain), and how many are left free (specifying the range). The 
domain is the tensor product of some number of Cj times some number of C* , and so 
is the range. In the extreme case that all incoming and outgoing chevrons on the spin 
network have kets and bras plugged into them the result is simply a number and the 
range is C. In that case the spin network, as a C-valued linear operator, can be seen 
as a vector in the space dual to the domain. 

This facile identification of closely associated operators, and their reinterpretation 
as elements of vector spaces, is an important advantage of spin networks. It is difficult 
and awkward to do something similar with the Dirac notation. 

In general, a tensor is a multilinear operator acting on a tensor product of some 
set of vector spaces and their duals. Thus, a spin network is a notation for a tensor on 
some product of the Cj and their duals. The edges of the spin network terminating in 
ket or bra chevrons indicate the nature of the space the tensor acts on. In general, the 
tensor product of two tensors in spin network notation is indicated by the placing of 
the two spin networks together on the page, in any orientation. The outer product of 
a bra and a ket illustrated in Fig. [7] is a special case. Partial or complete contractions 
of tensors are indicated by joining some or all ket chevrons with bra chevrons. 



2.5. The 2j symbol and intertwiner 

We now consider another intertwiner, which leads to an important mapping between 
kets and bras, alternative to Hermitian conjugation. The intertwiner acts on the 
Hilbert space Cj (g> Cj, the tensor product of two carrier spaces of the same j. In 
general we wish to consider the second carrier space as distinct from the first, which is 
the purpose of the prime on the second factor. To within a normalization and phase, 
there is a unique vector in this space that is invariant under rotations; we call it \K), 
and express it in terms of the Clebsch-Gordan coefficients by 

\K) = V27TT J2 \3 m ) ® \j m ') Cjjmm'- (1) 
mm' 

This vector can also be expressed in terms of the "2j-symbol," which we define in 
terms of the usual 3j-symbol by 

3 3 \_( 3 3 \ oo Mr"' 5 , (2) 

m m' J ~\m m' J ~ ~ ^/2]~+T m '~ m { ' 

The terminology "2j-symbol" is not entirely standard, but it has been used by 
Stedman (1990). The invariant vector \K) can also be written, 

\K) = J2 \3m) ® \jm') (-iy- m 5 m ^ m . = ]T \jm) ® \j, -m) (-l)^ m . (3) 

mm' m 

By Hermitian conjugation we convert \K) into the bra (K\, which is otherwise 
an intertwiner K : Cj ® ; —> C Just as the components of the intertwiner W are 
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Ill 



3 , J t 



m 



{-iy- m s ri 



m -^—%> 3 —^- J - <f^— m' 



Figure 10. Components of the map K, a standard bivalent intertwiner. 



the 3j-symbol, the components of the intertwiner K are the 2j-symbol, multiplied 
however by \/2] + 1 because of a normalization convention. Figure [10] shows first the 
spin network for the components of K, which is conceived of as a standard bivalent 
node or intertwiner. The small arrows indicate that the components on the first line 
can be considered the result of plugging basis kets into the intertwiner itself, as seen 
on the second line. The short line extending above the node is a "stub" (Stedman 
1990), whose purpose is to orient the node. The convention is that if we start at the 
stub and move in a positive (counterclockwise) direction, the first and second edges 
we encounter are respectively the first and second operands of K, conceived of as a 
map : Cj ® -> C. 

We note that if V and W are vector spaces, then V x W is not the same as V ® W, 
but if we have a bilinear map on V x W it can be extended to a linear map on V ® W 
by linear superposition. For example, in the previous paragraph we have regarded K 
as a bilinear map : Cj x Cj — > C, and computed its components as K(\jm), \jm'}). 
The first and second operands of this expression correspond to the first and second 
edges as specified by the stub. 

{K \ =) 3 i 3 < 



Figure 11. Spin network for map K : Cj ® — > C. 

Figure [10] also gives the numerical values of the components of K, and in the 
final diagram, the network for K itself, with two kets inserted. The network for K in 
isolation is illustrated in Fig. [Til regarded as a bra vector on Cj (g> Cj, that is, as an 
element oiC*®C*. 



m 3 I 3 . — m! = V2j + 1 



m 



in 



Figure 12. The stub can be regarded as a vestigial edge of a 3.7-symbol with 
value j = 0. 

The stub in Fig. [10] can be regarded as a vestigial edge of a 3j-symbol or W- 
intertwiner with the value j = 0, although one must beware of the normalization 
convention. This is illustrated in Fig. [T^] which is equivalent to ©. The value does 
not depend on the direction of the arrow on the zero edge (see below for rules for 
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reversing arrows). 

= (_i)2j ^, — , J i J ". — ^ = (_i)jy ^ — , J | J . — v 

Figure 13. The intertwiner A" acquires a phase (— l) 2 ^ if the two operands are 
swapped. 

The components of K : Cj® — > C, seen in Fig. [TU1 acquire a phase of (— l) 2 - 7 if 
m and m' are swapped. This is equivalent to the statement 

K{\4>)M)) = {-1)*K{\<V)M)), (4) 

for all \ip), \<f>), which is illustrated in spin network language in Fig. 1131 The final 
diagram differs from the preceding simply by a 180° rotation, so the value is the same. 
But this leads to the rule, that a spin network acquires a phase of (— l) 2 -' when the 
stub at a bivalent node is inverted. In particular, the arrow on the null edge in Fig. [T121 
can be inverted without changing the value. 

2.6. Kets to bras 

> — j i j — « j 4> = > — j \ j . v 

= * y- J —r 

Figure 14. The intertwiner K can be used to convert a ket to a bra. 

The intertwiner K : Cj ® Cj —> C has an alternative interpretation as a map 
K\ : Cj — > C* , where for simplicity we have dropped the prime on the second factor, 
and where the 1-subscript distinguishes the new map from the old. This alternative 
interpretation is natural in spin network language, as seen in Fig. 1141 The spin network 
for £ Cj is plugged into the second operand of the spin network for K, resulting 
in a spin network with one free bra chevron. The choice of the second operand of K 
for this purpose is conventional. The result is an element of C*, that is, it is a bra. 
As indicated in the figure, we abbreviate this spin network by drawing the same spin 
network for \tp) that we started with, except the ket chevron is converted into a bra 
chevron. 

In other words, when we use K\ to convert a ket into a bra, we just flip the ket 
chevron, leaving everything else the same. In particular, we do not put a star on 
the label of the ket. This distinguishes the map K\ : Cj — > C* from the metric or 
Hcrmitian conjugation, which is also a map : Cj — > C*. When the metric is used to 
convert a ket to a bra, not only is the ket chevron flipped to a bra chevron, but a star 
is appended to the label. These two maps are quite distinct; in particular, K\ is a 
linear map, while the metric is an antilinear map. As indicated in Fig. 1141 the results 
are not the same. 

When a ket is turned into a bra, the components with respect to some basis 
change from contravariant to covariant. But since there is more than one way to do 
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this, any notation based on the position (upper or lower) of the indices is inadequate 
to represent the result. 

4^- @ 

Figure 15. An arbitrary spin network with one edge terminating in a ket chevron. 

The map K\ can be used to turn a ket chevron into a bra chevron on any spin 
network, not only on kets themselves. A notation for an arbitrary spin network with 
an edge terminating in a ket chevron is shown in Fig. 1151 The circle around the X 
indicates the rest of the spin network, which may include other edges terminating in 
ket or bra chevrons. 

> — 3 . 3 — « 3 — Q = > — 3 . t Q 

= yJ. — Q 

Figure 16. Converting a ket chevron into a bra chevron with the map K\ : Cj — > 

q. 

To convert the ket chevron in Fig. [T5] into a bra chevron we simply insert it into 
the second operand of K : Cj <£> Cj — > C, as shown in Fig. [T()J The remainder of the 
spin network, indicated by the X, does not change. 

2. 7. Bras to kets 

< 3 1 3 > 

Figure 17. Spin network for the maps K~ x : C* ® C} -> C and K^ 1 : C* — > Cj. 

(j)* < (( 3 I 3 )) V* 

7 7 

= 0* — i i ^ — ^* 

Figure 18. Spin network for the quantity A" - 1 ((</>|, (i/'D- 

The map i^i : Cj — > C* has an inverse, the map K± 1 : C* — > C 3 , that takes bras 
into kets. We associate iff 1 with a closely related map K^ 1 : C* ® C* — > C that 
is expressed by the spin network in Fig. [TTJ This spin network and the meaning of 
the —1 on K^ 1 are to be defined, but the spin network represents a linear operator 
: C* ® C* — > C with the ordering of the two operands being specified by the stub. The 
action of K~ x on two bras is illustrated in Fig. [TBI 



Semiclassical Mechanics of the Wigner 6 j -Symbol 



12 



^ J ^ J , J y _ J , J y 

= 4,*— ^ ^_ 

Figure 19. Action of if^ -1 on a bra. 

We now define the network in Fig. [17] by requiring that K^ 1 act on a bra by 
inserting it into the first operand of that network, as illustrated in Fig. [TH and by 
requiring that K^ 1 actually be the inverse of K\. Figure fT9l shows the action of K± 1 
on an arbitrary bra (cf>\. The use of the first operand of K~ x for this purpose is a 
convention, but one that makes our overall notation for mapping kets to bras and 
vice versa consistent (see Fig. [22] below). As indicated, we abbreviate the result by 
taking the original network for the bra (4>\ and simply flipping the direction of the 
chevron. We do not unstar the identifying symbol. As indicated, the result differs 
from Hermitian conjugation applied to (<f>\, which is the ket \(f>) (without the star). 

< 1 I 3 » 3 1 3 < = < 3 < 



Figure 20. The requirement that K 1 o K = Id. 

The requirement that K^ 1 actually be the inverse of K\ is illustrated in Fig. l20l 
The stub on the network for K^ 1 is inverted so that the output of the first step is fed 
into the first operand of K~ x . 

3 j | 3 3 



Figure 21. Two arrow flips, with stubs oppositely oriented, annihilate each 
other. 

The identity represented by Fig. [201 is usually encountered in practice in the form 
shown in Fig. 1211 (a fragment of a spin network). The 2j-node inverts the direction of 
the arrow. Two such inversions, with stubs pointing in opposite directions, annihilate 
one another. 

^ I ; , J ^> J i J <^ J ' \ 3 ' y 

= ( 3 , I i i i I ' y = ^ 3 ' \ 3 ' y 



Figure 22. Flipping the chevrons on K gives K . 

We have made an independent definition of the spin network in Fig. 1171 but it is 
the same as the spin network for K, shown in Fig. 1 1 1 1 with both bra chevrons flipped. 
Since we now have a convention for flipping bra chevrons (by applying K^ 1 ), for 
consistency we must show that the two results are the same. This is done in Fig. 1221 
which uses the identity of Fig. [21] 
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m*— j - « 3 I 3 )) -^ m'* 

= m * 1 i I m >* = (-1)3-™ 6 m ,- m , 

Figure 23. The components of K -1 have the same numerical values as the 
components of K. 

By inserting resolutions of the identity into the diagram in Fig. [20] it is easy to 
work out the components of K~ x . These are displayed in Fig. [23] Notice that they 
have the same numerical values as the components of K (see Fig. ITO]) . 

,Q 3 « 3 . 3 > = Q — i \ j y 

Figure 24. Converting a bra chevron to a ket chevron on an arbitrary spin 
network. 

By using K~ x we can convert a bra chevron into a ket chevron on any spin 
network, not only on bras and kets themselves. This is illustrated in Fig. [24] which 
may be compared to Fig. [TBI 




Figure 25. Reversing the arrow on an edge of a spin network incurs a phase of 
(-1)*. 

Finally, by using Figs. HH [24] and [21] it may be shown that when we reverse the 
arrow on an edge of a spin network, we incur a phase of (— l) 2 ^ . This is done in Fig. [25] 

2.8. Raising and lowering indices 

When we convert a ket to a bra by the action of K\, then the bra has components 
with respect to the standard basis \jm) that are simple functions of the components 
of the original ket with respect to the standard basis (jm\. Mapping the one set of 
components to the other is "lowering the index." Using K± 1 to convert a bra to a 
ket similarly amounts to "raising the index." More generally the procedure can be 



Semiclassical Mechanics of the Wigner 6 j -Symbol 



14 



applied to an edge of any spin network terminating in a starred (contravariant) or 
an unstarred (contravariant) index. The index can refer to any basis, not just the 
standard one. 




Figure 26. Expressing contravariant components in terms of covariant 
components. 




Figure 27. Rules for raising and lowering indices. 

Figure [26] shows how to the express contravariant components in terms of the 
covariant components in the standard basis. By plugging in the numerical values of 
the components of K, we obtain the first line of Fig. [57J Similarly we derive the 
second line of Fig. [57] for expressing covariant components in terms of contravariant 
components. 



mi m * 




Figure 28. The completely covariant and completely contravariant components 
of the 3j-intertwiner are numerically equal. 



m — J I K — m! = m* ■ I I m'* = (-iy~ m 5 m _ m / 

Figure 29. The completely covariant and completely contravariant components 
of the 2j'-intertwiner are numerically equal. 

If these rules are used to raise all three covariant components of the 3j-intertwiner 
(Fig. [I), then we find that the completely contravariant components have the same 
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values, namely, the 3j-symbol. This is illustrated in Fig. [25J To show this it 
is necessary to use the symmetry of the 3j-symbol (see Varshalovich et al 1981, 
Eq. (8.2.4.6)). Then by setting one of the j's to zero and using Fig. [TU we find that 
the same is true for the completely covariant and completely contravariant components 
of the 2j-intertwiner, as shown in Fig. 1291 The same result is obtained by comparing 
Figs. M and [23] 

2.9. Hermitian conjugation of spin networks 

Consider a spin network of arbitrary complexity involving only 2j- and 3j-nodes. The 
network is allowed to have any number of edges terminating in bra or ket chevrons, or 
in starred or unstarred labels such as m indices. By using the identities above, possibly 
with the insertion or removal of 2j-nodes and the extraction of phases of the form 
(— 1) 3 , it is possible to bring the spin network into a standard form, in which all edges 
joining 3j-nodes have arrows pointing toward the 3j-node, all edges joining 2j-nodcs 
have arrows pointing away from the 2j-node, all edges terminating in a starred symbol 
have arrows pointing toward that symbol, and all edges terminating in an unstarred 
symbol have arrows pointing away from that symbol. Next, by inserting resolutions of 
the identity, which involve m-sums, it is possible to express the spin network as sum 
over the completely covariant components of 3j-symbols and completely contravariant 
components of 2^-symbols, times a tensor product of bras and kets. 

In this form it is easy to take the Hermitian conjugate. Under Hermitian 
conjugation, bras go to kets and vice versa, while the covariant components of 3j- 
symbols and contravariant components of 2j-symbols do not change, since they are 
real. By using Figs. [28] and [29] however, these components can be rewritten as the 
completely contravariant components of 3j-symbols and the completely covariant 
components of 2j-symbols. The m-sums can now be done, reversing the earlier 
insertions of resolutions of the identity. Then the other steps leading to the standard 
form can be reversed. 

The result is a simple rule for the Hermitian conjugation of any spin network 
of the given form: All ket chevrons are changed to bra chevrons and vice versa, the 
directions of all arrows are reversed, and all edges terminating in a symbol have a star 
added to the symbol, with a double star being removed. 

2.10. Discussion of spin network rules 

Our rules for spin networks differ from those of Yutsis et al and most of the literature 
in the Yutsis tradition primarily by our ability to express abstract vectors (kets), dual 
vectors (bras) and tensors in addition to the components of those objects. Also, we 
indicate the nature of an m index (covariant or contravariant) by the presence or 
absence of a star, rather than the direction of the arrow. One result is that our rules 
for reversing the direction of the arrow are more uniform than in the Yutsis tradition, 
where such a reversal picks up a phase (— l) 2 -? only on internal edges. In our approach, 
the rule applies everywhere, including edges terminating in an m index. In addition, 
our rules for Hermitian conjugation are simpler than those in the Yutsis tradition, 
where phase factors must be introduced. The simplification is due to the explicit 
introduction of 2j-symbols, and the use of stubs. 

To translate a Yutsis spin network into one of ours, it is necessary only to put 
stars on m indices terminating edges with outward pointing arrows. 
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The standard form of a spin network discussed in Sec. 12.91 has all 3j-nodes with 
inward pointing arrows, and all 2j-nodes with outward pointing ones. If we assume 
the standard form, then the arrows become superfluous and can be dropped. Only the 
stubs and 2j-nodes remain, in comparison to a Yutsis-style spin network. This is the 
procedure advocated by Stedman (1990). For the purposes of this paper we will keep 
the arrows, since we wish to have finer control on the spin network than that offered 
by the standard form. 

3. Models for the 6j-symbol 

For given values of the six j's, the 6j-symbol is just a number, but to study its 
semiclassical limit it is useful to write it as a scalar product (B\A) of wave functions 
in some Hilbert space. This can be done in many different ways, corresponding to what 
we call different "models" of the 6j-symbol. In this section we describe a class of such 
models that are related to one another. We begin by summarizing our notation for 
the Schwinger representation of angular momentum operators. Then we present what 
we call the "12j-model," which was used by Roberts (1999) in his derivation of the 
Ponzano-Regge formula. Next we describe the "4j-model" which we will use for the 
semiclassical analysis of this paper. We also mention an 8j-model for the 6j-symbol. 

3.1. The Schwinger Representation 

Our notation for the Schwinger representation of angular momentum operators is 
similar to that used in I. We denote the Schwinger Hilbert space by SS = L 2 (R 2 ); it 
is the space of wave functions ip( x i> x 2) for two harmonic oscillators of unit frequency 
and mass. The usual annihilation and creation operators are = (x^ + ip M )/-\/2, 
<^u = { x ij. ~ *i?/i)/v2, for fx = 1,2; we use hats on operators to distinguish them from 
their classical counterparts. We define operators 

/ = Ji = ^a)o- t a, (5) 

where i = 1,2,3 and where a (without the \i index) is a 2- vector (or column spinor) of 
operators, with a' the adjoint (or row spinor) and with obvious contractions against 
the Pauli matrices o~i. These operators satisfy the commutation relations [/, Ji] = 0, 
[Ji, Jj] — ieijk Jk- We set H = 1. Note that I = (Hi + H2 — l)/2, where H^, n = l,2, 

are the two harmonic oscillators. There is also the operator relation J = /(/ + 1). 
We denote the squares of 3- vectors in bold face. The operators J generate an SU(2) 
action on SS, which carries one copy of each irrep j = 0, 1/2, 1, . . ., that is, 

55 = ^©^-. (6) 

3 

The irreducible subspace Cj is an eigenspace of J with eigenvalue j. In the semiclassical 
analysis of spin networks, the spaces Cj that the spin networks refer to are interpreted 
as one of the irreducible subspaces of a Schwinger Hilbert space SS. Similarly, is 
interpreted as a subspace of a space SS* . In this way the bra and ket vectors referred 
to by the spin network are interpreted as wave functions on R 2 , and the spin network 
itself can be interpreted as a wave function in K 2Ar . 

In the various nj-models we take tensor products of the Schwinger Hilbert space, 
writing SS r for the r-th copy. Similarly, we put an r index on various operators, for 
example, d r , I r , J r , r = 1, . . . , n, or, with two indices, a rfi , ^ — 1, 2. 
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3.2. The \2j-model of the Qj-symbol 

We begin with the standard spin network (Yutsis et al 1962) of the 6j-symbol, shown in 
Fig. [501 According to the remarks in Sec. 12.101 this spin network can be reinterpreted 
according to our conventions, presented in Sec. [2] without modification. We will refer 
to the labeling of the six j's in the 6j-symbol shown in Fig. [30] as the "symmetric" 
labeling. 

h h is 1 / i \ 

k 35 36 J /jS^lK 

4 



Figure 30. The Yutsis spin network for the 6,j'-symbol, in a symmetrical labeling 
of the j's. The numbers 1, 2, etc on the spin network refer to ji, j2, etc. 

We perform an operation on each edge of the spin network that is illustrated for 
edge 1 in Fig. [3T] After the second equality edges are labeled by 1 and 1'. These 
refer to two distinct carrier spaces, say, Cj 1 and Cj x , with the same value of j (that 
is, ji). The introduction of such distinct carrier spaces does not change the value of 
the 6j-symbol, which is just a number. In the second equality we have expressed the 
lower ket chevron as a bra chevron transformed by K± , as in Fig. 1241 In the final 
diagram the arrows are directed toward both 3j-nodes connected by the original edge, 
and a 2j-node has been inserted. We do this on all six edges of the spin network in 
Fig. 1301 The resulting diagram is somewhat busy so we do not attempt to draw it, 
but each edge of the original diagram now looks like the final diagram in Fig. I3T1 




Figure 31. By inserting 2j-symbols in each edge of a Yutsis diagram, all 3j- 
symbols can be put into standard form (purely contravariant) . 



■ mi mj 



Ji 



Ji 



Ji 



Figure 32. Breaking lines to represent 6,j'-symbol as a sum over 2j- and 3j- 
symbols. 

We now break up the final diagram in Fig. [31] in two different ways. The first 
way is illustrated in Fig. 1321 in which the primed and unprimed lines are broken into 
summations over primed and unprimed quantum numbers m. This is done for all 
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six edges of the original spin network. The 6j-symbol is represented as a product of 



six copies of a 2j-symbol and four of a 3j'-symbol. 
2j-symbol, the result is 



Using the definition @ for the 
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m' R 



3i 

m.4 



J3 34 35 

m' 3 TO4 rri5 



(7) 



This formula may be compared to Edmonds (1960), eq. (6.2.3). Edmonds uses what 
he calls a "metric tensor" (really the components of K or if" -1 , multiplied by (— l) 2 - 7 , 
see his eq. (3.7.1)), which relative to our 2j-symbol introduces an overall phase of 
rir=i( — l) 23r - He also swaps m r and m' r for r = 4,5,6 relative to our definitions, 
which introduces further phases. The product of these phases is 1, showing that the 
formulas agree. 

Another way of breaking up the 6j-symbol is to stop with the third diagram of 
Fig. |3T] in the transformation of the six edges of the 6j'-symbol. Again the resulting 
diagram is too busy to draw, but it can be regarded as the complete contraction of two 
tensors, one the tensor product of six 2j-intertwiners, all terminating in ket chevrons, 
the other the tensor product of four 3j-intcrtwiners, all terminating in bra chevrons. 
There are twelve ket chevrons and twelve bra chevrons altogether, which we think of 
as living in twelve carrier spaces Cj r and C'^ and their duals, where r = 1, . . . , 6. These 
can be viewed as subspaces of twelve Schwinger Hilbert spaces, SS ri SS' r , r = 1, . . . , 6, 
and their duals. Then the 6j-symbol takes on the form (B\A), where states \A), \B) 
belong to the total Hilbert space Hi2j — (Y\ r ®SS r ) ® ([\ r ®SS' r ). These states are 
illustrated in Fig. ESI where the state \B) has been turned into ket form by Hermitian 
conjugation. 



\B) = 



1' 



2' 



3' 



1' 



\A) = 



A\ A\ A\ A\ A\ 



3' 



4' 



SK SI/ N/ M' M/ 



6' 



Figure 33. The two states that appear in Roberts' (1999) expression for the 
6j'-symbol. 

The usual custom in physics is to specify a state by the operators and quantum 
numbers of which the state is a simultaneous eigenstate. This requires that the 
eigenstate be nondegenerate, so that it is determined to within a normalization and 
phase. This in turn requires, in a certain sense, that the number of independent 
operators should be equal to the number of degrees of freedom of the system. We 
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will not attempt to be precise about this statement, but will illustrate the principle 
in several examples. 

One example is the first 3j-state appearing in Fig. 1331 which lies in the Hilbert 
space SS\ <g> SS2 <8> SS3 and is a simultaneous eigenstate of I r , r = 1,2,3 with 
eigenvalues j r , r = 1, 2, 3. It is also an eigenstate of the vector of operators 

J 123 =J 1 +J 2 +J 3 , (8) 

the total angular momentum on this Hilbert space, with eigenvalue 0. That this 
simultaneous eigenstate is nondegenerate follows from standard angular momentum 
theory; and the number of operators (six) equals the number of degrees of freedom in 
the Hilbert space (two for each SS r , r = 1,2, 3). Thus we write this state as illustrated 
in Fig. [Ml indicating both operators and eigenvalues. This is otherwise the state \ W), 
illustrated in bra form in Fig. [B] As for the normalization and phase, these must be 
supplied by context. For the 3j-state illustrated in Fig. [34j these are given in terms 
of the 3j-symbol by Fig. [T] in particular, the state is normalized. 



h h I3 J123 
h h h 




Figure 34. Ket notation for invariant state of the 3j-typc. 



Similarly, the first 2j'-state in Fig. [33J lies in the Hilbert space SSi <E> SS[ and is a 
simultaneous eigenstate of the operators I\ and I[ with eigenvalues ji and ji , as well 
as of the total angular momentum operator on this space, 

3ii'=Ji+Jl, (9) 

with eigenvalue 0. This state is a simultaneous eigenstate of five operators, but in a 
sense only two of the three components of 3i±i are independent, so we should count 
only four independent operators, which agrees with the number of degrees of freedom 
in the Hilbert space (again, two each for SS\ and SS[). We write this state as 
illustrated in Fig. 1351 the normalization and phase are given by the components of 
K^ 1 shown in Fig. [23l In particular, with the square root factor in Fig. [35l this state 
is normalized. 

/|\ 

V 

h i'x Jii' \ = 1 _ 

h h / v/2ji + 1 

1 



Figure 35. Ket notation for invariant state of the 2j'-type. 



In this notation we can write the equations of Fig. [33J in the form 



\B) 



I\ I2 h J 123 

h h 33 



I'l I5 I§ Jl'5'6 
31 35 36 



^2 ^6 I* J2'6'4 
32 35 36 



I3 I4 h J3'4'5 
33 34 35 



(10) 
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and 

\a) = ( n v2>+t 

\r=l 

One can see that Edmonds' form of the 6j-symbol ([7]) is equal to (-B|A). 

This scalar product was the starting point for Roberts' (1999) analysis of the 
asymptotics of the 6j-symbol. We shall comment below on further aspects of Roberts' 
calculation. 



h A J u , 
3i h 



h I'e, J66' 



(11) 



5.5. XTie Triangle and Polygon Inequalities 

We make a remark on a generalization of the triangle inequalities before presenting 
the 4j-model of the 6j-symbol. If (^1,^2,^3) are three nonnegative lengths, the usual 
triangle inequalities are \ii — £j\ < £k < li + lj, where (i,j,k) = (1,2,3) and cyclic 
permutations. We generalize these as follows. Let {li, i = 1, . . . , n} be a set of lengths, 
li > 0, i = 1, . . . , n. Then this set is said to satisfy the "polygon inequality" if 

1 " 

max{4}< ( 12 ) 

i=l 

This is equivalent to the triangle inequalities when n = 3. In general, it represents the 
necessary and sufficient condition that line segments of the given, nonnegative lengths 
can be fitted together to form a polygon with n sides (in R N , N > 0). 



3.4. The Aj -model of the 6 j -symbol 

A different way of writing the 6j-symbol as a scalar product begins with Fig. 1361 in 
which the j's in the 6j-symbol of Fig. [30] have been relabeled. We will refer to the 
labeling in Fig. [36] as the "asymmetric" labeling, which is more appropriate for the 
T/-model. After the relabeling, we have reversed the arrow on the edge j'3 of the spin 
network (labeled simply by 3), incurring a phase (— l) 2 -* 3 , and broken four edges into 
scalar products of a bra and a ket. Then, on the second line of Fig. [36] we have 
unfolded the bra and the ket, written the bra as the Hermitian conjugate of a ket, and 
adjusted phases. The result expresses the 6j-symbol as a phase times a scalar product 
of two states lying in the Hilbert space — Yl r= x ®SS r . It is understood that the 
ket terminating a line labeled r lies in the j r -irreducible subspace of SS r . 

These same states arise in the recoupling of four angular momenta with a resultant 
of zero. Let Z be the subspace of Cj 1 ® Cj 2 <£> Cj 3 ® Cj 4 upon which 
4 

Jtot=^Jr = 0. (13) 
r=l 

This is the subspace of rotational invariants, that is, Z is the set of states \ip) € Cj 1 <E) 
Cj 2 (g) Cj a C3> Cj 4 that are invariant under rotations. According to the rules for addition 
of angular momenta, subspace Z is nontrivial (dim Z > 0) if J2t=i 3r = integer and if 
the set {j r } satisfies the polygon inequality (fT2j) . In accordance with the remarks in 
Sec. 12. 3[ the subspace Z can also be interpreted as the space of 4-valent intertwiners, 
that is, S'J7(2)-invariant maps : Cj ± ® Cj 2 (g) Cj 3 ® Cj 4 — > C. The notation Z is a 
mnemonic for "zero" (the eigenvalue of Jtot)- 

As far as recoupling theory is concerned the spaces Cj r can be any carrier spaces 
of SU(2) for the given values of j r , but in our application we shall interpret the space 
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Figure 36. Asymmetric labeling of the 6j'-symbol, and decomposition of spin 
network. 

Cj r as the irreducible subspace j r of the r-th copy of the Schwinger Hilbert space SS r . 
Then Z becomes a subspace of H^j. We shall assume that the fixed values of the 
four j r , r — 1, . . . , 4 are chosen such that dimZ > 0. We should properly label Z by 
the four j r values since T-L^j contains as many subspaces of the type Z as there are 
choices of the four j's. For simplicity, however, we will suppress this dependence in 
the notation, it being understood that j r , r = 1, . . . , 4 are given. 

Standard recoupling theory gives three ways of constructing an orthonormal basis 
in Z. One uses Clebsch-Gordan coefficients to couple angular momenta according to 
the pattern 1 + 2 = 12, 12 + 3 — 123, 123 + 4 = 0, resulting in the normalized state 
\B) lying in Z, expressed in terms of a spin network in Fig. [37] A second way couples 
according to the pattern 2 + 3 = 23, 1 + 23 = 123, 123 + 4 = 0, producing the 
normalized state \A) illustrated in Fig. [33 A third way, which we will not consider 
further, uses the intermediate coupling 1 + 3 = 13. 



\A) = (-l)h-h-h+h ^ 2 j 12 + 1 



B) = (-l)2(ii+J2+i4)^2j 23 + l 




Figure 37. Two ways of recoupling four angular momenta with a total of zero. 

The quantum numbers j\2 and j'23 of the intermediate angular momenta range in 
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integer steps between the bounds 

.712, min < jl2 < J 12, max, 
.723, min < J23 < .723, max, 

where the maximum and minimum values are given in terms of the four fixed j r , 
r = l,...,4by 

Jl2,min = max(|ji - j 2 |, \ja - H\), Jl2,max = min(ji + j 2 ,.73 + ji), ^ 
J23,min = max(|j 2 ~ h\, \jl ~ J*4 1 ) j J23,max = min(j 2 + js, jl + J4). 

Then the dimension of Z is given by 

D = dim Z = jl2, max ~ jl2,mm + 1 = .723, max ~ .723, min + 1. (16) 

An expression for dim Z can be given that is symmetrical in (ji, j 2 , J3, ji)- Using 
the fact that |x| = max(a;, —x), the difference between ji 2ima x and ji2,min becomes the 
shortest distance between one set of four numbers, {ji — j 2 , j 2 — ji, J3 — 3Ai3a ~ J3}i 
and another set of two numbers, {ji + j 2 , % + ji]- But this is the minimum of the 
distance between all eight possible pairs taken from the two sets. Thus 

D = dimZ = 2min(ji, j 2 , j 3 , H, s - ji, s - j 2 , s - j 3 , s - j 4 ) + 1, (17) 

where s is the semiperimeter, 

s = \{ji+32+h+h)- (18) 

More precisely, if D computed by (fT7| is < 0, then subspace Z is trivial (dimZ = 0); 
otherwise D — dimZ. The formula (fT7|) bears an interesting relationship to the Regge 
symmetries of the 6j-symbol (Varshalovich et al 1981, Eq. (9.4.2.4)). 

Now by comparing Figs. l36l and 1371 we see that the 6j-symbol is proportional to 
a scalar product, 

(B\A) = V(2ii2 + l)(2j 2 3 + 1) { {} f f j- (19) 

t J3 J4 J23 J 

This is the scalar product that we shall use for the semiclassical analysis of the 6j- 
symbol in this paper. In a different notation, we can write (J23 IJ12) instead of (B\A), 
which emphasizes the fact that this is a unitary matrix element connecting two bases 
on Z. The usual orthonormality relations satisfied by the 6j-symbol (see Edmonds 
eq. (6.2.9)) are equivalent to the unitarity of (J23 1 J12) - 

Notice that in this 4j-model, the angular momenta J r , r — 1,2,3,4 are 
independent operators, while the two remaining angular momenta, 

Jl2=Jl+J 2 , J23=J 2 +J 3 , (20) 

are not, rather they are functions of the first four. As usual in the Schwinger 
representation, the quantum number j r , r = 1, . . . , 4 specifies the eigenvalues of both 
I r and J^, that is, j r and j r (jr + 1), respectively. And the quantum numbers ji 2 
and j 2 3 specify the eigenvalues of the operators 3\ 2 an d J 2 3j that is, ji 2 (ji 2 + 1) and 
J23O23 + 1), respectively. But there are no operators Ji 2 or 7 2 3. 

The states \ A) and \B) in Fig. [37] can be expressed as eigenstates of complete sets 
of operators, 

h h h I a J 12 JtotA I m 



1^4) 



Jl ]2j3jA Jl2 



1 1 h h h J 23 Jtot 

jl 32 J3 34 323 , 



(21) 
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in a notation like that used in (|10p and As mentioned, these states are normalized, 
and their phases are specified by Fig. |3"T1 

Notice that each state \A) and |B) has a list of eight independent operators 
(counting the three components of Jtot), corresponding to the eight degrees of freedom 
in the 4j-model. We will call these lists of operators the A- list and -B-list, and write 
them collectively as 

A = (h,l2,h,h, J12, Jtot), . 

B = (Ii, I2, hi lit J23' Jtot)) 

We denote elements of these lists with subscripts, for example, Ai or Bi, i = 1, . . . , 8. 
The operators in the either of the lists (|2"2"|) do not commute with one another (because 
the components of Jtot do not commute; otherwise all commutators are zero), but they 
do possess simultaneous eigenstates in Z, which are unique to within a phase, namely, 
the states (|2Tj). 

3.5. The 8 j -model of the 6 j -symbol 

We obtain an 8j-model of the 6j-symbol by inserting 2j-symbols into edges 12 and 
23 of the spin network of Fig. [35] and then treating them in the same way as the 
2j-symbols in the 12j-model. The result is the Hilbert space Hsj = (Ilr=i ®SS r ) (g> 
SSi2(8>SS' 12 <& SS23 ( S)SS' 23 . This model has more symmetry than the 4j-model but less 
than the 12j-model. Operators 7 12 and I23 exist in this model (as well as operators 
I[ 2 and I23), unlike the 4j-model. This gives the 8j-model certain advantages over 
the 4j-model. We shall not consider the 8j-model further in this paper. 

4. The Classical Manifolds 

In this section we study the classical mechanics that will be relevant for the 
semiclassical analysis of the 6j-symbol in a 4j-model. We begin by presenting our 
notation for the Schwinger phase space and products of it that are used to represent 
coupled, classical angular momenta. Other spaces that will be important are obtained 
by Poisson and symplectic reduction. Then we examine the geometry of the A- 
and B-manifolds in phase space that support the states (|2"TT) . including a rather 
general analysis of why they are Lagrangian. Finally we discuss the Bohr-Sommcrfcld 
quantization of these manifolds. 

4-1. The Schwinger Phase Space and Other Spaces 

The classical phase space for two harmonic oscillators (the Schwinger phase space) is 
$ = (M 4 ,<ip A dx), with coordinates (x,p) € K 4 , for x,p £ K 2 . See Sec. 1.3 for more 
details on our use of the Schwinger representation, as well as Cushman and Bates 
(1997) for a more detailed discussion of the geometry of two harmonic oscillators and 
the role played by the Hopf map. Here dp A dx means Y] , t — 1 dp^ A dx^ 7 and similarly 
for other obvious contractions over fj, = 1, 2. Complex coordinates = (x^ + ip^) / 'V%, 
Zfi = (x^ — ip^/^/2 are also useful, allowing us to write $ = (C 2 ,idz A dz). 
Interesting functions on $ (classical observables) include 

I = -zz, Ji = \zuiZ, (23) 
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where i = 1,2,3, obviously the classical analogs of ([3]), and where z or z without 
indices indicates a 2-component "spinor," that is, an element of C 2 . In comparison 
to © notice the absence of the hats, indicating that these are classical observables. 
These satisfy the Poisson bracket relations {I,Ji} = and {Ji,Jj} = Cijk Jk, as well 
as the identity J 2 = I 2 . The Hamiltonian flow of / is a £7(1) action on $, while 
the flows of Ji, i = 1,2,3, generate an SU(2) action on $. Both actions are easily 
expressed in the complex coordinates z: that of U(l) is z H > exp(—itp/2)z, where ip is 
the variable conjugate to I, while that of SU(2) is z n- uz for u € SU(2). 

The orbit of the U(l) action generated by I passing through any point z^Oon 
$ is a circle on which tp is a coordinate, covered once when < ip < 4ir. The level set 
/ = J for J > 0, where I : $ — > R is the function and J 6 R is the contour value, is a 
3-sphere to which the orbits of I are confined. The definition (|2"31 of Jj is interpreted 
as a map ir :<!>—> R 3 , where the three components J$ of J are coordinates on R 3 (thus, 
this space is "angular momentum space" ) . This map is a Poisson map (Marsden and 
Ratiu 1999), giving R 3 the Poisson structure { Ji, Jj} = e^k Jk- We denote M 3 with 
this Poisson structure by A, as indicated by the diagram, 

(J = J) C $ 

7T]i TT (24) 

E c A 

The map tt can also be interpreted as the momentum map (Abraham and Marsden 
1978) of the SU(2) action on $, so that A or angular momentum space is identified 
with su(2)*. When tt is restricted to a level set I = J > in $ (a 3-sphere), it projects 
onto the 2-sphere |J| = J in A. This is the projection map tth of the Hopf fibration, in 
which the orbits of / are the fibers (or "Hopf circles"). It is also the map of symplectic 
reduction (Abraham and Marsden 1978) of the level set / = J in $ by the U(l) action, 
so that the 2-spheres |J| = J in A arc symplectic manifolds. The symplectic form on 
one of these 2-spheres is J dfl, where dtt is the element of solid angle (to within a 
sign). We denote one of these spheres with its symplectic structure (for some value 
J > 0) by E; these spheres are also the symplectic leaves of the Poisson structure in 
A. 

In an nj-model of a spin network we take Cartesian products of the Schwinger 
phase space $ to obtain the phase space for n independent classical angular momenta. 
We will illustrate the notation for the 4j'-model. We write = $ 4 for the entire 
phase space; apart from the symplectic structure, this is C 8 = R 16 . Coordinates on 
$4j are x^ r , p^ r , z^ r , etc, fj, = 1, 2, r = 1, . . . , 4. We denote the r-th copy of $ by 
$ r , and define functions I r , J r , on <I> r on the pattern of (|23p . that is, just by adding r 
subscripts to all the variables in those equations. Naturally these can also be viewed 
as functions on $4^ . The vector 3 r is the r-th classical angular momentum. 
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We generalize the diagram (|24l) to the 4j-model as follows, 

(I r = Jr) C $ 4 j 



n 



CL 



C 



c 



A 



(25) 



r 



where tt means dividing $4^ by the £7(1) 4 action generated by I r , r = 1, . . . , 4. Thus 
J r , r = 1, . . . , 4 are coordinates on A4J = (R 3 ) 4 , in which the Poisson bracket of two 
functions / and g is given by 



Similarly, for four positive contour values J r > 0, r = 1,...,4, the level set 
I r = |J r | = ,/,., r = 1,...,4, in $4^ is (S 3 ) 4 (indicated simply by I r — J r in the 
diagram). The map ir restricted to this space is a power of the Hopf map (simply 
denoted tth in the diagram), which projects (S* 3 ) 4 onto the space £4., — (S 2 ) 4 , a 
symplectic manifold in A4J . Notice that the radii of the spheres (the 3-spheres in $4., 
and the 2-spheres in Aij) need not be equal for different r. The spaces CL and T in 
(I25|) will be explained later. 

Other important classical observables on $4j are J12 = Ji + J2 and J23 = J2 + J3, 
their squares, 3\ 2 and J23, and the total angular momentum J to t = J2r=i ^r- 

4-2. Level Sets and Contour Values 

The A- and i3-lists of operators in (|22|) correspond to lists of classical observables 
(without the hats), 



We will denote the members of these lists by Ai, Bi, i = 1, . . . , 8. As discussed in 
I and in Littlcjohn (1990), the Lagrangian manifolds that support the semiclassical 
approximations to the states \A) and \B) are the level sets in $4^ of these lists of 
classical observables, with quantized values of the contour values. We will defer the 
question of quantization to Sec. 14.61 and for now just examine these level sets for some 
suitable contour values. We will call the level sets of the A- and £?-lists the A- and 
£?-manifolds. 

It will be convenient to distinguish notationally the functions in the two lists, 
regarded as maps : $4^ — > IR, from the contour values, which are real numbers. Our 
notation is summarized in Table [TJ The vector of functions Jtot is given the value 
because that is the only contour value we will consider. The other contour values 
are variable. The conventions in the table solve some notational problems in I, but 
the notation requires care. For example, the magnitude of the vector J r , regarded 




(26) 



A — {h,h, h, -^4, Ji2: Jtot), 
B = {h,l2,h,Ii, J23, Jtot)- 



(27) 
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function 


value 


It 




T 2 

J 12 


7 2 

■'12 


T 2 

J 23 


7 2 

J 23 


Jtot 






Table 1. Notation for functions : <J>4j — > R and values (real numbers) of those 
functions. In the first row, r = 1, . . . , 4. 

as a function : $4^ — > R, is not J r (because J r is a number, not a function), but on 
the level set I r = J r it is true that that |J r | = J r , in view of the identity between 
functions, J 2 = Z 2 . If we wish to refer to the magnitude or the square of the vector 
J r , regarded as a function on $, we will write |J, | or J 2 (in bold face), not J r or J 2 . 
We will denote the lists of contour values collectively by a or b, so that 

a = (Ji, J2, J3, Ji, J12, 0), 
6 = (J l5 J 2 , J3, J4, J 23 , 0), 

with components <Zj, 6j, i = 1, . . . , 8. 

Consider now the conditions on the contour values a, b such that the A- 
and B-manifolds should exist (as nonempty sets of points). The question of the 
dimensionality of these manifolds will be postponed to Sec. 14.41 but it turns out that 
their maximum (for any contour values) and generic dimensionality is 8, one half the 
dimension of the 4j Schwinger phase space. We work with the A-manifold, since the 
conditions for the B-manifold are completely analogous. 

First, the A-manifold clearly does not exist unless J r > 0, r = 1, . . . , 4, so let us 
assume this condition. Then if the ^-manifold exists, we can pick a point on it and 
evaluate the vector functions J n r = 1, . . . , 4, to determine the projected point in A4J. 
We visualize this point as four vectors in a single copy of R 3 . In addition we compute 
the vector J 12 = Ji + J 2 and plot it along with the others in R 3 . Next, we move 
the vectors (Ji, J2, — J12) end-to-end to create a triangle with sides (Ji, J2, Ji2)- We 
move the vectors by parallel transport in R 3 , that is, without rotating them. Thus, 
the triangle inequalities are satisfied in the triplet of lengths (Ji, J2, ^12)- Also, since 
Jtot = 0, we can form a second triangle out of vectors (J3, J4,Ji2), so the triangle 
inequalities are satisfied in the triplet (J3, J4, J12). The two triangles have the edge 
J12 in common, creating a figure like Fig. 1381 

Conversely, suppose that J r > 0, r = 1, ... ,4 and that the triangle inequalities 
are satisfied in triplets {J\, J2, J12) and (J3, J4, J12). This means that two triangles 
can be constructed in IR 3 with the given lengths. By translating and/or rotating the 
two triangles, we can bring the 12 edges into coincidence, thereby creating a figure 
like Fig. [38l and hence a point of A^j. But since the map tt : $4^ — » is onto, there 
exists an inverse image of this point in $4^, hence the A-manifold exists. 

In summary, we have shown that the A-manifold exists iff J r > 0, r = 1, . . . , 4, 
and the triangle inequalities are satisfied in the triplets ( J\, J2, J12) and (J3, J4, Ji2)- 
A similar result holds for the B-manifold, where the triplets are ( J2, J3, J23) and 
(Ji, J4, J23)- The triangle inequalities imply that J12 and J23 lie in the ranges, 

</l2,min < Jl2 < ^12, max, (29) 
^23 ,min < J23 < ^23 (30) 




Jl2 



Figure 38. A point of the A- or £?-manifolds in <E>4j = C projects under 7r onto 
a point of A47 = [M 3 ] 4 on which the four angular momenta J r add to zero. We 
visualize this point as four angular momentum vectors in a single copy of R 3 that 
add to zero. In such diagrams, the vectors need not be based at the origin; here 
we place them end-to-end. Also shown in the figure is the vector J12 = Ji + J2- 

where 

Ji2,min = max(| Ji - J 2 \, I J3 - Ji\), Ji2,max = min(Ji + J 2 , J3 + Ji), (31) 

^23, min = max (I J 2 - J 3 |, \Ji - Ji\), J23, max = min(J 2 + J 3 , J\ + Ji). (32) 

These inequalities imply J12, J23 > 0. Notice that the value is not always excluded, 
for example, if J\ = J 2 and J3 = J4 then Ji 2 = is allowed. 

It turns out that the A- and B-manifolds, when 8-dimensional, are Lagrangian. 
According to the Liouville- Arnold theorem (Arnold 1989) the compact level sets of 
complete sets of Poisson commuting observables are generically Lagrangian tori. The 
sets of observables A or B of interest in this paper are not commuting, however, so the 
Liouville- Arnold theorem does not apply. In a similar situation in I we showed that 
the manifold in question (what we called the "Wigner manifold") was nevertheless 
Lagrangian. In the following section we present a general set of circumstances in 
which Lagrangian manifolds are obtained, which cover not only the cases considered 
in the Liouville-Arnold theorem but also all cases we know of in the asymptotics of 
spin networks, including the manifolds studied in I and the A- and B-manifolds of 
this paper. 

4-3. Level Sets, Orbits, and Lagrangian Manifolds 

We use general notation in this section that differs somewhat from that of the 
application to the 6j-symbol in the rest of the paper. The basic conclusion of this 
section is the following. Let {Ai, i = 1, . . . , m} be a collection of classical observables 
on a phase space (symplectic manifold) P of dimension 2N, that is, Ai : P — > M. 
Let the set {Ai} form a Lie algebra, that is, the Poisson brackets {Ai, Aj} are linear 
combinations of the Ai. Let L be the level set Ai = ai, for some contour values ai, 
and suppose that all Poisson brackets {Aj, Aj} vanish on L. Finally, suppose L is a 
smooth manifold of dimension (thus, m> N). Then L is Lagrangian. The reader 
who is willing to accept this conclusion can skip the remainder of this section. 

We deal with classical Hamiltonian systems with symmetry. Basic references 
on this subject are Abraham and Marsden (1978), Marsden and Ratiu (1999) and 
Cushman and Bates (1997). We make quite a few assumptions in this section, but 
most of them are generic. The most important one that is not is the assumption that 
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the momentum, in the general sense of the value of the momentum map, is a fixed 
point of the coadjoint action of the group, that is, it is G-invariant. 

Let P be a symplectic manifold of dimension 2N, let G be a connected Lie group 
of dimension m with Lie algebra g, dual Q*, and symplectic action on P. Suppose 
the momentum map M : P — > g* exists, that is, the action of G is generated by 
Hamiltonian flows of a set of Hamiltonian functions. Let {£i,i = 1, . . . , m} be a basis 
in g, and let c k ij be the structure constants, so that = c k ij£k (summation 

convention). Define functions Ai : P —> R by Ai(x) = (M(x),£i) for all x £ P. These 
functions form a Lie algebra under the Poisson bracket, {Ai, Aj} = c k ij Ak, with the 
same structure constants as the £j in g. Let Xj = cu~ 1 dAi be the Hamiltonian vector 
fields associated with the Ai, where u> is the symplectic form on P, regarded as a map 
from vector fields to 1-forms. 

For some point Xo G P let a = M(xq), so the level set L of the A's passing 
through xq is given by A = a (that is, L = M~ (ct)). For simplicity we assume 
that L has only one connected component, or else we restrict consideration to the 
connected component passing through x$. Let n be the rank of the set of differential 
forms {dAi,i = l,...,m}, assumed to be constant over L. Then n < m and 
dimi = 2N — n. Also let the orbit of the G-action through xq be B. The vectors 
{Xi, i = 1, . . . , m} are tangent to B and span its tangent space at each point. Also, 
rank{A"i} = rank{cL4i} = n, since u) is nonsingular. Thus, dim_B = n. 

At this point we have the basic geometry for symplectic reduction, as illustrated 
in Fig. [39] Shown in the figure is the intersection I of L and B, which is the orbit of 
xq under the isotropy subgroup of a under the coadjoint action of the group. Dividing 
L by the isotropy subgroup produces the reduced symplectic manifold. 





/b 


I 

/ 






L / 









Figure 39. Schematic illustration of the basic geometry of symplectic reduction. 
L is the level set of the momentum map passing through xo, B is the orbit of the 
group action, and I the intersection, also the orbit of the isotropy subgroup. 

Now suppose that a £ g* (the generalized momentum) is a fixed point of the 
coadjoint action of the group, that is, Ad*a = a, Vg G G. By differentiation and 
contraction with an arbitrary element of g this implies (a, [rj, £]) — 0, Vr/, £ € g, or, 
by setting T) = £j, £ = £j, c k ijOk — 0. But this implies {Ai, Aj} — c k ij Ak — on 
the level set L, where A^ — ak- In other words, the functions Aj, which may form 
a nontrivial Lie algebra on P, have vanishing Poisson brackets among themselves on 
L. This in turn implies that the Ai's are constant along each other's flows on B, that 
is, XiAj = — {Aj, Ai} = 0, so B C L. Therefore dimS < dimL, or n < 2N — n, or 
n < N. 
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Finally, let us assume that n takes its maximum value n = N. Then dim B = 
dimi = N. If L is compact, then B = L, that is, the level set and the orbit 
coincide. Moreover L is Lagrangian, since the vectors Xi span the tangent space 
to L and uj(Xi,Xj) = — {A;, A,} = on L. This in turn implies that L supports 
locally a solution of the simultaneous Hamilton- Jacobi equations involved in finding 
a semiclassical eigenfunction of the operators with principal symbols Ai. 

In the case of an integrable system the group is Abelian, G — M. , the N classical 
observables Ai Poisson commute everywhere in phase space, and on a generic level set 
A = a the differentials dAi are linearly independent everywhere. Since the group is 
Abelian, a is automatically a fixed point of the coadjoint action, and thus the level 
set is Lagrangian. Moreover, if it is compact, it is an iV-torus. 

In the case of the manifolds explored in I (the jm- and Wigner manifolds) and 
the A- and _B-manifolds of this paper, the isotropy subgroup of the action of G on 
xq E P is zero-dimensional (a discrete set of points), so dimG = dimB = m = n. 
Moreover, n = N, so the manifolds are Lagrangian. 

4-4- Properties of the A- and B -manifolds 

We now apply the analysis of Sec. l4.3l to the Amanifold. In the following we visualize 
a point of $4j as a 2 x 4 matrix of complex coordinates z^ r , referring to the r-th 
column as the r-th "spinor." 

Here the symplectic manifold is P — $4j, and we take the group to be 
G = U{lf x SU(2), where the action of U(l) 5 is generated by (h, I 2 , 1 3 , h,3l 2 ), 
and that of SU{2) is generated by Jtot- For given r, the observable I r generates a 
{7(1) action, the multiplication the r-th spinor by e _l ^ r / 2 , where ip r is the evolution 
variable conjugate to I r , with period 4-7T. See (1. 27). Also, Jtot generates the diagonal 
SU(2) action, that is, all four spinors are multiplied by the same element u £ SU(2) 
(the 2x4 matrix is multiplied on the left by u). See (1.28) and (1.29). As for the 
observable Jf 2 , Hamilton's equations are 



where t is the variable of evolution. See (1.28) for a similar calculation, and notice 
that 3 12 is a constant of the flow generated by 3\ 2 . Thus, Jf 2 generates a U(V) action 
that rotates spinors 1 and 2 about the axis ji2 = J12/IJ12I, that is, it multiplies them 
by u(ji2,9) G SU(2) (in axis-angle notation for an element of SU(2), see (1.18)) while 
leaving spinors 3 and 4 invariant. The period is 8 = 4tt or t = 2tt/J 12 . The U(l) 
action generated by 3\ 2 is not a rotation about a fixed axis, since J12 depends on 
where we are in $4^, but it does commute with the other £7(1) actions generated by 
the I r , as well as the SU(2) action generated by J. Similar statements can be made 
about J 23. 

Now we determine the maximum dimensionality of the ^4-manifold. Since there 
are 8 functions in the A-list and dim $4^ = 16, the answer will be 8 if the functions 
are independent. Functions are independent at a point if their differentials, in this 
case, {dAi,i = 1,...,8}, are linearly independent. In any case, the rank of this 
set of differentials is the same as the rank of the set of Hamiltonian vector fields 
Xi = Lo~ 1 dA i . 

The A-manifold cannot have its maximum dimension if any J r — 0, for r 
1 1 . since that implies dl r — 0. So we assume J r > 0, i = 1, . . . , 4. This means 




ifr = l,2, 
if r = 3,4, 



(33) 
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that all four spinors z r are nonzero, and hence the orbit of the group f/(l) 4 generated 
by the I r , r = 1, ... ,4 is a four-torus T 4 . This is also the fiber of the projection 
7r : $4j -t A 4j . 

The condition J r > 0, r = 1, . . . , 4 also means that all four vectors J r , r = 1, . . . , 4 
are nonzero. We think of these vectors in a single copy of R 3 . When the group SU(2), 
whose action is generated by Jtot, acts on $4^ , the effect on the vectors 3 r , r = 1, . . . , 4, 
is to rotate them by the corresponding element of 50(3) (see (1.21)). If any of these 
vectors is moved by the 50(3) rotation, then in $4^ we will have moved off the initial 
fiber of the projection tt, hence we will have motions that are linearly independent 
of the U(1) A action. To make all three independent directions of rotation in 50(3) 
give rise to linearly independent motions, we require that the isotropy subgroup of the 
50(3) action on the set {J r , r — 1,...,4} (that is, the diagonal action) be trivial. 
This is at a point at which J to t = Yl r =i = ^- This requires that the four vectors 
J r , r — 1, . . . , 4 (which are nonzero) be noncollinear (otherwise the isotropy subgroup 
is SO(2) or the whole group 50(3)). This in turn requires that at least one of the 
triangles, ( J\, J2, J12) or (J3, J4, J12), have nonzero area. 

Next there is the [/(l)-action on $4j, generated by 3\ 2 . This has the effect on 
the vectors J r , r = 1, . . . , 4 of a rotation of vectors Ji and J2 about axis J12, while 
leaving J3 and J4 invariant. In order that this motion move us off the initial fiber, 
we require that Ji and J2 be linearly independent, that is, that triangle ( Ji, J2, J12) 
have nonzero area. And in order that this motion be linearly independent of overall 
rotation generated by Jtot, we require that the triangle (J3, J4, J12) also have nonzero 
area, in order that the rotation of Ji, J2 about the axis J12 should change the shape 
of the figure created by the four vectors. (The overall rotations generated by Jtot do 
not change the shape of the figure.) 

We conclude that at a point of the A- manifold, rank{cL4.i, i = 1, . . . , 8} = 8 iff 
J r > 0, r = 1, . . . , 4 and none of the triangle inequalities in ( Ji, J2, J12) or ( J3, J4, J12), 
or equivalently, the inequalities in (|29]l . is saturated. Notice in particular that this 
implies J12 > 0. Since these conditions depend only on the contour values and not 
where we are on the A-manifold, the A-manifold, under the stated conditions, is a 
smooth, 8-dimensional manifold. 

In addition, the A-manifold is compact. It is also connected, as follows by 
consideration of its projection under tt (see Sec. 14. This means that the action 
of the group U (l) 5 x SU(2) is transitive on the A-manifold (it is the orbit of any point 
on it under the group action). It then follows from the discussion of Sec. 14.31 that 
the A-manifold, when 8-dimcnsional, is Lagrangian. Similar statements apply to the 
B-manifold. 

To find the topology of the A-manifold, we first find the isotropy subgroup of the 
action of U(l) 5 x SU(2) on a point xq on this manifold. If we denote coordinates 
on U(l) 5 x SU(2) by (0i, ip2, 1P3, V>4, 0, u), where u e SU(2) and where the five 
angles are the 47r-periodic evolution variables corresponding to (Ji,J2,/3,/4,J? 2 ), 
respectively, then the isotropy subgroup is generated by two elements, say, x — 
(2-7T, 2tt, 27r, 2tt, 0, —1) and y — (0, 0, 2n, 2ir, 2n 7 — 1). The isotropy subgroup itself 
is the Viergruppe {e, x, y, xy} = (Z 2 ) 2 . Thus the A-manifold is topological^ 
U(l) 5 x SU (2) /(Z 2 ) 2 . This is the same logic used in I to find the topology of the 
"Wigner manifold" of that paper. The analysis is the same for the i?-manifold, which 
has the same topology. The isotropy subgroup would be larger in degenerate cases, 
for example, when some triangle inequalities are saturated. 

Now it is easy to find the invariant measure on the A- or B-manifolds. It 
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is dipi A dip2 A ^"03 A dip4, A d9 A du, where du is the Haar measure on SU(2) 
(du = sin f3 da A df3 A dj in Euler angles). Thus, the volume of the A- or B- manifold 
with respect to this measure is 



where the 1/4 compensates for the 4-element isotropy subgroup. 
4-5. Projections and tetrahedra 

We now study the projection of the yl-manifold onto A^-, for contour values such 
that the manifold is 8-dimensional. The projection consists of the set of four nonzero 
vectors in R 3 , {J r , r = 1, . . . , 4} of given lengths J r > with a sum of zero, Jtot = 0, 
creating a closed chain of links as in Fig. [38j such that J12 has a given, positive value 
and both triangles 1-2-12 and 3-4-12 have nonzero area. We choose to place the four 
vectors end-to-end in the order (Ji, J2, J3, J4), as in Fig. 1381 this is an arbitrary choice, 
but convenient for studying observables 3\ 2 an d J23 (if we wished to examine J^ 3 we 
would choose a different order). By filling in the lines J12 and J23, the closed chain 
of links becomes a tetrahedron, as shown in Fig. |4"01 



The set of all such figures, modulo orientation, is the circle 5 1 — 50(2), whose 
angle is the dihedral angle about the edge J 12. Including the orientation, we see that 
the projection of the A-manifold onto A 4 j is topologically 5*0(2) x 5*0(3). One might 
imagine that this should be described as an 50(3) bundle over 50(2), but in fact 
the bundle is trivial. This is connected, so the inverse image under 7r, which is the 
A-manifold itself, is a 4-torus bundle over 50(2) x 50(3), and is also connected. An 
action of the group 50(2) x 50(3) on A^y is generated by 3\ 2 and J to t, regarded 
as observables on A±j using the Poisson structure (|26|) . This is the projection of the 
action of U(l) 5 x SU(2) on $4j. The motion generated by 3\ 2 is the "butterfly" 
motion illustrated in Fig. |40l in which the butterfly flaps one of its wings. 

As the triangle 1-2-12 rotates around the edge 12, the length | J 2 3 1 varies between 
a maximum and minimum value. The minimum value is reached when triangles 1-2-12 
and 3-4-12 lie in the same plane on the same side of line 12, and maximum when on 
opposite sides of line 12; these are also the minimum and maximum of | J23 1 on the 
A-manifold in $4^. These extremal values are reached when the tetrahedron is flat 




(34) 




Figure 40. The function J^ 2 generates a rotation about axis J12 of vectors Ji 
and J2, rotating triangle 1-2-12 by an angle 9 while leaving triangle 3-4-12 fixed. 
The result is a family of tetrahedra with values of j'23, the length of vector J23, 
ranging from a minimum to a maximum. 
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(its volume is zero), and can be obtained in terms of the other five J's as the roots of 
the Cayley-Menger determinant (Berger 1987, Crippen and Havel 1988), 
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where x — | J 2 3 1 . Such determinants were an important part of the analysis of Ponzano 
and Regge (1968). In this case the determinant expresses the condition that the volume 
of the tetrahedron is zero. Expanding the determinant gives a quadratic equation in 
x 2 in terms of the five J's specifying the A-manifold. The same condition can be 
expressed in terms of a smaller (3 x 3) Gram matrix of dot products, as discussed in 
Littlejohn and Yu (2009). 

A similar analysis applies to the £?-manifold except that here we fix the length 
J23 of vector J23, and create tetrahedra of different shapes by varying the dihedral 
angle between triangles 4-1-23 and 2-3-23, that is, by rotating vectors J2 and J3 about 
the axis J23 = J23/J23- In this process, the length IJ12I varies from some minimum to 
some maximum, which can be obtained by replacing x 2 in (f35|) by Jf 3 , and then J 2 2 
by x 2 where now x — | J 12 1 - 



4-6. Quantizing the Manifolds 

The A- and i?-manifolds can be subjected to Bohr-Sommcrfeld quantization, which 
selects certain contour values as quantized. The process uses the Weyl symbols or 
transforms (Weyl 1927, Wigner 1932, Groenewold 1946, Moyal 1949, Berry 1977, 
Balazs and Jennings 1984, Littlejohn 1986, Ozorio de Almeida 1998) of the lists (|22|) 
of operators, which are functions on $4^. The Weyl symbols of selected operators 
are summarized in Table [2] It is important that the classical manifolds be the level 
sets of the Weyl symbols of the operators whose simultaneous eigenfunctions we seek 
(Littlejohn 1990); this holds in the present case because the Weyl symbols in the table 
are always equal to the corresponding classical observable (without the hat), to within 
an additive constant. Notationally we could have defined the classical observable 
(without the hat) as the Weyl symbol of the corresponding quantum observable, but 
the conventions in the table make it easier to establish the connections with the usual 
quantum numbers in physics. In particular, the zero point energy has been subtracted 
from the quantum observables I r , but not from the classical ones I r , which explains 
the 1/2 on the first row of the table. 

To quantize the A- or B-manifold, we first determine the homotopy group, then 
we compute action integrals and Maslov indices along generators of the group, then 
we require that the action plus Maslov correction be an integer multiple of 2tt. Only 
manifolds of full dimensionality (8) can be quantized. This is the procedure followed 
in I, and the analysis is very similar in this case; in particular, in both cases the 
homotopy group is Abelian. 

We just summarize the results, speaking of the ^4-manifold. We find that J r , 
r = 1, . . . , 4 are quantized in half- integer steps, which we write in terms of the quantum 
numbers j r as 

J r = j r + 1/2, (36) 
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J12 - 3/4 
J 2 3 - 3/4 



Table 2. Weyl symbols of selected operators. In rows containing operators I r , 
J r and J%, r = 1, . ..,4. 

where the allowed values of j r are 0, 5, 1, • • •• Smaller values of j r are not allowed 
because for j r = —1/2 the manifolds do not have full dimensionality, while for 
j r < —1/2 they do not exist. We choose the quantum number j r so that it agrees 
with the usual notation in physics for the eigenvalues of various operators, but that is 
not confirmed until we compute the semiclassical eigenvalues in Sec. 14.71 

Similarly, J12 must be an integer or half-integer on quantized manifolds, which 
we write in terms of a conventional quantum number by J12 = j'12 + 4. In addition, 
there is the condition that j\ + J2 +J12 and j'3 + j'4 + J12 be integers, the requirement 
(with the usual interpretation of the quantum numbers) that the 3j-symbols in Fig. [36] 
should exist. See also eq. (1.85). These imply that j% + ji + J3 + j'4 must be an integer, 
part of the conditions that the subspace Z defined in Sec. 13.41 be nontrivial. 

The range of the quantum number j'12 is determined by the requirement that the 
A-manifold be 8-dimensional. Looking first at the upper limit, if J12 is quantized we 
have 

•hi = jn + \< J 12, max = min(ji +j 2 + l,j 3 +J4 + !)■ (37) 

Given the other conditions on j'12, this implies that the maximum quantized value 
of J12 is Ji2.max — \- Similarly, we find that the minimum quantized value of J12 
is Ji2,min + \- The quantized values of J12 are separated from the maximum and 
minimum classical values (for given J r , r = 1, . . . , 4) by a margin of ^, and are spaced 
in integer steps. These rules imply 

Jl2,max = >/l2,max _ 1; Jl2,min = ^12,min- (38) 

They also imply the bounds (|15[) on the quantum number jyi- Similar results apply 
to the B-manifold and the quantized values of J23. 

Figure [4T1 is a numerical example of these quantization rules. The square in the 
figure is given by the bounds ([29)) and (|30|) . while the spots are the quantized values 
of J12 and J23. The latter are related to the usual quantum numbers by 

1 1 

Jl2=jl2+2> J23=j23 + -. (39) 

The spots form a square array because O23 IJ12} is a square matrix. The size of the 
matrix (the number of rows or columns of spots) is dim given by (1161) . Other 
features of this figure will be explained later. 



Semiclassical Mechanics of the Wigner 6 j- Symbol 34 

Jl2 




Figure 41. The J12-J23 plane for ( Ji, J2, J3, J4) = (5,^,6,-^) (all quantized 
values). Region U is where either the A- or B- manifold does not exist. Inside the 
square both manifolds exist; A and F are the classically allowed and classically 
forbidden regions, respectively. Caustic curve is shown. 

Note that the minimum value of J r , r = 1, . . . , 4 and of Ji 2 on any quantized A- 
manifold is 5 , so the corresponding vectors always have a positive length. In addition, 
the triangles 1-2-12 and 3-4-12 always have a positive area. Similar conclusions apply 
to the quantized B-manifolds. 

4-7. Semiclassical Eigenvalues 

Once the classical manifolds are quantized, we find the semiclassical approximations to 
the eigenvalues of the operators in the A- or B- lists by evaluating the Weyl symbols of 
those operators on the classical manifold. Doing this for the operators I r , r = 1, . . . , 4, 
we find the eigenvalue of I r is j r ; this is the exact answer, something that was to be 
expected since I r is a quadratic polynomial in the fundamental x's and p's of the 
system, and Weyl quantization of such operators is exact. See Littlejohn (1986) 
for more on classical and quantum quadratic polynomials. If we now compute the 
eigenvalue of J 2 by using the operator identity J 2 = I r (I r + 1), naturally we get 
the exact answer, j r (j r + !)• On the other hand, if we evaluate the Weyl symbol 
of the operator J 2 on the quantized level set, we get (according to the table) the 
eigenvalue (j r + 1/2) 2 — 3/8 = j r (jr + 1) — 1/8, with an error of —1/8. There is an 
error because J 2 is a quartic polynomial in the fundamental x's and p's of the system, 
so Weyl quantization is not exact. The error is, however, of relative order h 2 , that is, 
1/j 2 , which is the error expected in lowest order semiclassical approximations. The 
operators J 2 , do not appear in the A- or B-lists, but operators J 2 2 and J 2 ^ do, and 
here again there is an error at order . Moreover, in this case there are no operators 
J12 or I23 which could be used to obtain the exact eigenvalues. This is a drawback 
of the 4j-model in comparison to the 8j- or 12j-models, where such operators exist. 
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It does not, however, change any of the subsequent analysis, which depends only on 
using the quantized manifolds to carry out the stationary phase calculation. 



4-8. Time Reversal 

The action of the antilinear time reversal operator on a carrier space Cj can be 
defined by 

e\jm) = (-iy- m \j,-m) (40) 

(Messiah 1966). This is equivalent to = Ky 1 o G, where G is the antilinear metric 
or map of Hermitian conjugation, and K^ 1 is defined in Sec. 12.71 The composition 
of the antilinear G with the linear K^ 1 is the antilinear time reversal map 0. The 
map is easily extended to the full Schwinger Hilbert space SS and tensor products 
thereof such as HAj ■ 

Classically the antilinear : Hij — ► Hij corresponds to an antisymplectic map 
: $4j — > $4j. In the complex coordinates, its action on all four spinors is given by 

( Zlr ) = ( ~~ Z2r ) , (41) 



%2r J \ z l 

that is, : z r 1— > vz r , where z r G C 2 and where 

v = = (° 1 ~ 1 V (42) 

We recall that in quantum mechanics, it is time reversal, not parity, that reverses 
the direction of angular momenta. At the classical level, this means that J r (0(x)) = 
—3 r (x), where x £ <&4j. 

Time reversal can be projected via tt onto A4J , where its effect on the coordinates 
is : J r i- > 3 r . It is an anti-Poisson map on A^j. 



5. Intersections and Actions 



In this section we consider the intersections of the A- and _B-manifolds, assuming 
that J r > 0, r = 1,...,4 are given. For now we treat this as a classical problem in 
which the J's (including J12 and J23) are continuous variables, but we note that if 
J r , r — 1, ... ,4 are quantized then they are automatically positive. The motivation, 
however, is to find the stationary phase points of the scalar product (flU)) . which are 
the intersections of the quantized manifolds. 



5.1. Classically allowed and forbidden regions 

A simple analogy will help to understand the results. Consider a one-dimensional 
harmonic oscillator, H = (1/2)(.t 2 +p 2 ) (classical or quantum). For a given value of 
the energy E, the classically allowed region is the interval of the x-axis between the 
turning points, given by x — ±y/2E, while the classically forbidden region is outside 
this interval. The classically allowed region can also be defined as the region of the 
x-axis where the two curves in phase space, x — const and H = E, have intersections. 
These curves are level sets of the observables x and H, which appear on the two sides 
of the matrix element when we write the energy eigenfunction as iPe{x) — (x\E). 
Inside the classically allowed region the intersections between the two curves consist 
of two points, with opposite momentum values. These are related by time reversal 
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(p — > —p). In the classically forbidden region the two curves have no real intersections, 
but if we complexify phase space and the curves x = const and H = E (maintaining 
real contour values x and E), then they do have complex intersections that are related 
to the exponentially decaying wave function in the classically forbidden region. 

We can view the classically allowed and forbidden regions in the x-E plane, in 
which both x and E are variables. See Fig. [42] The darkly shaded region in the figure, 
labeled UU, is the region E < 0, for which the level set H = E does not exist. Above 
the line E = 0, both level sets, x = const and H = E exist. The unshaded region, 
labeled AA, is the classically allowed region, while the lightly shaded region, labeled 
FF, is the classically forbidden region. The quantized values of the energy (E = n+^) 
are indicated as spots on the E-axis. There is only a one-dimensional array of spots 
because the observable x has a continuous spectrum. The parabola E = x 2 /2 is the 
caustic curve, separating the classically allowed from the classically forbidden region. 



5.2. Intersections of the A- and B -manifolds 

Similarly, for positive J r , r = 1, ... ,4, either the A- or B-manifold does not exist if 
J 12 or J23 lies outside the bounds ([25)) and ([ST)]). Those inequalities define a square 
region of the J23-J12 plane, as illustrated in Fig. [UJ A figure like this was first given 
by Neville (1971). The darkly shaded region outside the square, labeled U, is where 
either the A- or B-manifold does not exist. In the interior of the square both manifolds 
exist and have full dimensionality. On the boundary of the square they exist, but have 
less than full dimensionality, since some triangle inequality is saturated. If the A- and 
B- manifolds exist and intersect, then we are in the classically allowed (unshaded) 
region in the interior of the square. If they exist but do not intersect, then we are in 
the classically forbidden (lightly shaded) region. The caustic curve is the oval curve 
in the figure, separating the classically allowed from the classically forbidden regions. 
Notice that it touches the square boundary at four points. 

If we fix a value of J12, we can regard the classically allowed and forbidden regions 



E 




x 



Figure 42. Regions of the x-E plane for the harmonic oscillator: UU where 
curves H = E do not exist, and AA and FF, the classically allowed and forbidden 
regions, respectively. The parabola is the caustic curve. 
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as intervals of the J23-axis. The interval within which the i?-manifold exists is given by 
(|3U]k this is the interval of allowed J 2 3 values, given J r , r = 1, ...,4. The classically 
allowed region, on the other hand, is the interval of allowed J23 values, given J r , 
r = 1, . . . , 4 and J12 (equivalent to the statement that the manifolds intersect). Since 
this is a more restrictive condition, the classically allowed region must be a subset of 
the interval ([30]) . Moreover, since the A-manifold is connected, the subset must be a 
connected interval, since this subset is the range of J23 values that occur on a given 
^4-manifold. This is just what we see in Fig. |4TJ for most values of J12, the interval 
inside the square consists of a classically allowed region, surrounded on both sides by 
classically forbidden regions, outside of which are the regions in which manifolds of 
the given J23 values do not exist. The bounds of the classically allowed region are 
given by the roots of that is, they correspond to tetrahedra of zero volume. The 
caustic curve in Fig. @J]is the contour of detM = 0, where M is the Cayley-Menger 
determinant (pJ5j) . 

The volume vanishes if the tetrahedron is flat. For most values of the parameters, 
this does not require that any of the triangular faces have zero area. But for certain 
values of J12, the volume vanishes at a place where the face 2-3-23 has zero area (thus, 
J2, J3 and J23 are linearly dependent). This is the point where the caustic curve 
touches the boundary of the square on the left or right. Similarly, the caustic curve 
touches the boundary of the square on the top and bottom of the square, where not 
only does the volume vanish, but also the area of triangle 1-2-12. 

The caustic curve never lies outside the bounds ([30]) defined by the triangle 
inequalities. Biedenharn and Louck (1981) appear to claim the contrary, but there 
is the question of whether one is talking about the classical or quantum triangle 
inequalities. That is, the caustic curve does pass outside the bounds given by the 
square array of quantized spots. It also seems to us that the interpretation of 
Biedenharn and Louck of Fig. 6 from Ponzano and Regge is incorrect. 

If the values of J r , r = 1, ... ,4 are quantized, then we can plot the quantized 
values of J12 and J23 as a square array of spots, as in the figure. See other 
comments on this array in Sec. 14.61 In Littlejohn and Yu (2009) we incorrectly stated 
that the quantized values of J12 and J23 can fall exactly on a caustic, citing the 
theory of Brahmagupta quadrilaterals (Sastry 2002). That theory shows that plane 
quadrilaterals with integer sides and integer diagonals exist, that is, flat tetrahedra 
with all integer edges. However, to represent a 6j-symbol, the sums of the integers 
around the faces of triangles must be odd, and this condition cannot be met. (The 
integers divided by 2 are the values of J r — j r + |, and the sum of the j r around 
the faces must be an integer.) A correct proof of the nonexistence of flat, quantized 
tetrahedra, credited to Adler, is given in a brief citation by Ponzano and Regge (1968). 

We now examine the intersections of the A- and i?-manifolds in greater detail. 
First, the A- and i?-manifolds (assumed to exist) have an intersection in $4, iff their 
projections onto A^ intersect. Furthermore, the intersection in $4^ is the lift of 
the intersection of the projections in A^, with a T 4 fiber over every point. These 
statements follow from the fact that over every point x € A^ there is a 4-torus fiber 
in $4j, and that if x lies on the intersection of the projections of the A- and B- 
manifolds, then the 4-torus belongs to both the A- and B-manifolds in $4^. This is 
the same logic used in I under similar circumstances. 

To find the intersections of the projections in A^j we require four vectors J r , 
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r = 1 , . . . , 4 that satisfy 



4 



| J r | — 



(43) 



|Ji+J 2 |= J: 



|Ja + J3 1 — ^23i 



for the given values of J r > 0, r = 1,...,4 and of J12 and J23. A nice way of 
constructing these vectors is given in Appendix A of Littlejohn and Yu (2009), which 
uses the singular value decomposition of the Gram matrix of dot products associated 
with the Cayley-Menger determinant (|35p . This method not only gives an explicit 
solution for these vectors at any point in the classically allowed region, it also shows 
that they are unique to within the overall action of 0(3). It is obvious in any case 
that if J n r = 1, ... ,4 is a solution of (l4U|) . then so is S3 r for any S £ 0(3). This 
method was generalized to the 9j-symbol in Haggard and Littlejohn (2010). 

The group 0(3) is conveniently decomposed into proper rotations in SO(3) and 
spatial inversion, which is time reversal in the present case. It is a basic fact of 
geometrical figures in R 3 that spatial inversion is not equivalent to any proper rotation 
unless the figure is planar. Another fact is that the orbit of a geometrical figure 
under 50(3) is diffcomorphic to 5*0(3) itself, unless the dimension of the figure 
is < 1. These issues are discussed by Littlejohn and Reinsch (1995, 1997) in the 
context of molecular configurations. They imply that except at the caustics, where 
the tetrahedron is flat, two tetrahedra related by time reversal are not related by 
any proper rotation. Therefore, except at the caustics, the solution set of (|43l) in 
A4J consists of two disconnected subsets, each diffeomorphic to SO(3), related by 
time reversal. Each subset consists of tetrahedra of nonzero volume related by proper 
rotations. At a generic point of the caustic curve, where the tetrahedron is flat but still 
2-dimensional, the two subsets merge into one, which is still diffeomorphic to SO(3). 

The intersections in $4., are the lifts of these intersections in A^y. Therefore, 
except at the caustics, the intersection of the A- and B-manifolds consists of two 
disconnected subsets, related by time reversal, where each subset is a T 4 -bundlc 
over SO(3). These subsets are 7-dimensional, so the A- and i?-manifolds, which 
are 8-dimensional, intersect in two 7-dimensional submanifolds. The situation can be 
visualized as in Fig. [43J where A n B = 1\ U 1% and where 1\ and I2 are the connected 
intersection sets, related by (see ([4"Tjl ). Each intersection set is an the orbit of the 
group U(l) 4 x SU(2), where U(l) 4 represents the phases of the four spinors and SU(2) 
is the diagonal action (thus, the group is generated by I r , r — 1, . . . , 4 and Jtot)- 




I2 



Figure 43. A schematic illustration of the intersection of the two 8-dimensional 
A- and B-manifolds in <E>4j . The two intersections I\ and I2 are 7-dimensional, 
and are related by time reversal. Paths for computing action functions Sa and 
Sb relative to initial points ao and bo on the two manifolds are shown. 
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The isotropy subgroup of this group is Z 2 , generated by element 
(2-7T, 27r, 2-7T, 27r, — 1), in coordinates (V>i, '02, i>3, V>4, u ) for group T 4 x SU(2), where 
u G SU(2). The volume of intersection manifold I\ or / 2 with respect to the measure 
dipi A dip2 A d?/>3 A di/ji A du is 

V> = ^(4tt) 4 x 16tt 2 = 2 n 7r 6 . (44) 

An interesting aspect of the method of Appendix A of Littlejohn and Yu (2009) 
for finding the vectors J r , r = 1, . . . , 4 is that it also works in the classically forbidden 
region, where it produces complex 3- vectors that satisfy (|4"3"]l . These solutions are 
determined modulo the action of spatial inversion and 50(3,C). In fact, as discussed 
in that reference, the vectors can be chosen so that two components are real and one 
purely imaginary, so that the symmetry group of the solution set is the Lorentz group 
SO{2, 1) C 5*0(3,0). Roberts (1999) and others have referred to the tetrahedra in 
the classically forbidden region as living in R 3 with a Minkowski metric, while those 
in the classically allowed region live in R 3 with a Euclidean metric. This is a correct 
interpretation of the situation for the 6j'-symbol, but we do not think it is appropriate 
for generalizations to other spin networks. For example, in the 9j-symbol (Haggard 
and Littlejohn 2010) the vectors in the classically forbidden region cannot be chosen 
so that two components are real, instead they belong to C 3 , and the symmetry group 
of the solution set is SO(3, C), not some Lorentz subgroup thereof. This simply means 
that to explore the complex Lagrangian manifolds, complex Euler angles must be used 
when following the Hamiltonian flows generated by Jtot • The complexified A- and B- 
manifolds in the 6j-symbol support the asymptotic forms given by Ponzano and Regge 
(1968) in the classically forbidden region. 



5.3. Actions and phases on the A- and B -manifolds 

Let the x-space wave functions associated with states \A) and \B) of ([21]) be denoted 
4>a{x) = (x\A) and i/)b(x) = (x\B), where x € R 8 . The semiclassical approximations 
to these wave functions involve phases e lSA ^ and e lSB ^ x \ where actions Sa(x) and 
Sb (x) are integrals of p dx = Ylru Pru dx r ^ from some initial points on the two 
manifolds to some final point. The initial points on the A- and i?-manifolds are 
denoted ao and bo, respectively, in Fig. 1431 they determine the overall phases of the 
states | A) and \B). 

As explained in I, the branches of the stationary phase evaluation of the matrix 
clement (B\A) are associated with the intersections of the A- and £?-manifolds, in 
this case the manifolds I\ and I2 , so the asymptotic form of the 6j-symbol has two 
branches. Moreover, the phase associated with each branch is Sa — Sb evaluated 
on the corresponding intersection manifold (see, for example, (1.3) or (1.5)), and is 
independent of where we evaluate it on that manifold. Figure |43] illustrates two points 
Pi and P2, on intersection manifolds I\ and I2, respectively, with paths that may be 
used for computing the actions Sa and Sb- 

In the following we shall be interested in the relative phase between the two 
branches. Let us define 

Si = Sai — Sbi, S2 = Sa2 — Sb2, (45) 

and 

rP2 rP2 rPi rPi r 

S = S2-S 1 =l -/ -/ +/ pdx=j>pdx, (46) 

J ap J hi J ao J bo J 
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where the final integral is taken along the path that goes from pi to P2 along the A- 
manifold and then back to pi along the £?-manifold. The relative phase is independent 
of the initial points ciq and 6 , and is moreover a symplectic invariant. The relative 
phase is easier to determine than the absolute phases of either branch, which are 
related to the overall phase convention for the 6j'-symbol. In computing the relative 
phase, we note that the loop integral in (|46|) can be evaluated with respect to any 
symplectic 1-form, such as the complex one used in I (see (1.64)), 

f /;P r v dx r ^ — Im j> 2J z rfJ _ dz Tl j,- (47) 

The loop integral can be transformed by Stokes' theorem into an integral of the 
symplectic form over the enclosed area, since on $4, = C 8 all cycles are boundaries. 

5. 4- Closing the loop in A±j 

We shall construct the closed loop giving the relative phase between the branches 
according to (|4"o]) by following the Hamiltonian flows of various observables. Let us 
define the signed volume of the tetrahedron as V = (1/6) Ji • (J2 x J3), and let us 
take manifold I\ to be the one on which V < 0, so that V > on 12- Time reversal 
changes the sign of the volume when mapping I\ into 1% and vice versa. Let us start 
at a point p of 1%, as in Fig. [44] Then by following the J^-flow we trace out a path 
that takes us along the A-manifold to a point q of I2 . We cannot use any of the other 
seven observables defining the A- manifold for this purpose, namely, (Ji, 1%, I3, J4, J), 
since their flows confine us to the intersection manifold 2^. 




Figure 44. The J^ 2 -flow takes us from a point p of intersection manifold I± along 
the A-manifold to point q of manifold I2', and then the j| 3 -now takes us back to 
point p' of Ii . To close the loop it is necessary to connect p 1 to p inside I\ . 

We see that the J^ 2 -flow actually does take us to manifold I2 by considering the 
projection of the path p-q in Fig. 23] onto A^ . We visualize the projected path as a 
transformation applied to a set of four vectors in R 3 that define a tetrahedron. The 
situation is illustrated in Fig. US In part (a) of that figure, we have four vectors J r 
that sum to zero, defining a tetrahedron of negative volume. The lengths J r > 0, 
r = 1, ... ,4, J12 and J23 are assumed to have the prescribed values, and vector 3 12 
is drawn (but not J23). We take the point p of Fig. 031 to lie on the T 4 fiber above 
this tetrahedron. The J^-flow rotates the 1-2-12 triangle about the 12-axis by the 
right-hand rule while leaving the 3-4-12 triangle fixed (see (|33j) ). that is, it rotates the 
1-2-12 triangle into the foreground. Let the angle of rotation be 20i2, where </>i2 is the 
interior dihedral angle of the tetrahedron along edge 12 in its original configuration. 
This brings triangle 1-2-12 through triangle 3-4-12 to the opposite side, creating a new 
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tetrahedron with the same lengths (the new J23 is the same as the old one) , hence the 
same dihedral angles, but with the opposite signed volume. The result is illustrated 
in part (b) ol Fig. 03] a tetrahedron that is the projection of a point q E I 2 in $ 4j , as 
illustrated in Fig. Thus we see that the J^-flow does take us from l\ to I2 along 
the A-manifold, as claimed. 



Having reached point q G I2, we can go back to I\ along the B- manifold 
by following the J^-flow, reaching point p 1 € 1\ as illustrated in Fig. The 
transformation in is illustrated in Fig. [45] Part (c) of that figure is the same 
as part (b), except that vector J23 is drawn and J12 is suppressed. The j| 3 -flow 
rotates triangle 2-3-23 about the axis j23 = J23/J23, while leaving triangle 1-4-23 
fixed. Let the angle of rotation be twice the interior dihedral angle along edge 23, 
that is, 2023, a s illustrated in part (c) of Fig. [35] The result is part (d) of that figure, 
a tetrahedron in which the volume has been inverted a second time, taking us back to 
the original (negative) volume in part (a). We arrive at point p' £ Ii, as in Fig. [44] 

It is clear that p' is not the same as the original point p, because if it were, the 
orientation of the tetrahedron in part (d) of Fig. [45] would be the same as that in part 
(a). Thus to create a closed loop in $4^, we must follow some path in 1\ taking us 
from p' to p, as in Fig. 0H 

We create this path from p' to p in 1\ in two steps. First we apply an ST/ (2) 
transformation to all four spinors at point p 1 , that is, a diagonal transformation, 
whose projection onto A±j is an S*0(3) transformation of the tetrahedron in part (d) 
of Fig. [45] returning it to the original orientation in part (a). This is a proper rotation 
of all four vectors J r , that is, a rigid rotation of the entire tetrahedron, and it is the 
final step in a cycle of rotations that transform the original tetrahedron in part (a) 





Figure 45. A cycle of rotational transformations that takes a tetrahedron in 
back into itself. 
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into itself. That figures (a) and (d) must be related by some proper rotation is clear, 
since the lengths of the sides are the same and the signed volume is the same. In 
fact, it is easy to see that the axis of the final rotation is J4, since this vector is left 
invariant by both the Jf 2 - an d J^-flows, and is the same in parts (a) and (d). As it 
turns out, the final rotation has axis — J4 = — J4/J4 (notice the minus sign) and angle 
204, twice the internal dihedral angle along edge 4. This is illustrated in part (e) of 
Fig. 05l which is the same as part (d) except that all vectors are drawn. The effect of 
the final rotation is illustrated in part (f), which is the same as part (a) except that 
all vectors are drawn. 



5.5. Angle of the final rotation 

To obtain the angle of the final rotation about axis — J4, we use the fact that the 
product of two reflections is rotation. Let a reflection about a plane P be Q{P). 
There are four ways to draw the angle between the planes, one of which is denoted by 
a in Fig. 1461 For a given choice of dihedral angle a, let the outward pointing normals 
of the two planes be n and m, as in the figure. Then 

Q(n)Q(m) = R(a,2a), (48) 

where we use axis-angle notation for the rotation R and where the axis a of the rotation 
is along the line of intersection of the two planes and is given by 

n x m n x m . 

|n x m| sma 

Although the angle a, the normals m and n, and the axis a depend on which of the 
four choices is made for the dihedral angle, the resulting rotation does not (although 
it does depend on the order in which the reflections are applied). 

11 



m 



Figure 46. The product of two reflections is a rotation. 

Let us denote the first rotation taking us from part (a) to part (b) of Fig. [45] 
in axis-angle form by i?i2(ji2, 2(/>i2), where the 12-subscript on R indicates that this 
rotation only affects vectors Ji and J2. It is obvious from Fig.|45]that the effect of this 
rotation on all the angular momentum vectors is to reflect them in the plane 3-4-12, 
that is, 

j; = Q(3-4-12)J r , (50) 

where the prime refers to the values of the vectors after the first rotation and where 
r = 1, . . . , 4. This applies for r = 1, 2 because the reflection in the 3-4-12 plane has 
the same effect as the rotation, and for r — 3,4 because the reflection does nothing 
to these vectors and the rotation does not apply to them. Note that the 3-4-12 plane 
is the same in parts (a) and (b) of the figure (the plane is not affected by the first 
rotation) . 
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Similarly, we denote the second rotation, taking us from part (c) to part (d) of 
Fig. SSI by R 23 (j23> 2( h3), where the prime on J23 indicates that the axis is the 23- 
direction after the first rotation. Then the effect of the second rotation on all four 3' r 
is the same as a reflection in the plane l'-4-2'3, that is, 

J.;' = Q(l'-4-2'3)j;. (51) 

The planes of the two rotations intersect in edge 4 of the tetrahedron, which can be 
seen more clearly in Fig. \$7\ which is the same as part (c) of Fig. 05] except that all 
vectors are drawn. In Fig. [47] plane 3-4-12 is the back plane, and is also the plane of 
the first reflection. Plane l'-4-2'3 (primes are omitted in the figure) is the plane of the 
second reflection. In comparison to (j48| . if we identify a with the interior dihedral 
angle ^4 at edge 4, then m is the outward normal to plane 3-4-12, while n is the 
outward normal to plane l'-4-2'3. Their cross product is in the direction J4, so we 
have 

J' r ' = i?(j 4 ,20 4 )Jr, r = l,...,4. (52) 

This is the rotation taking us from part (a) to part (e) of Fig. US] to undo that rotation, 
we apply R(— J4, 2</> 4 ) to pass from part (e) to part (f) of that figure. 




Figure 47. The plane of the first reflection is the back plane, 3-4-12, while that 
of the second is the plane l'-4-2'3 (primes suppressed in the figure). 

To summarize the rotational history, we have applied the rotations 

R(-U, 204)i?23(j23, 20 23 )i?12(jl2, 2012) (53) 

to the tetrahedron in part (a) of Fig. 05] taking it through a cycle of tetrahedra and 
returning it to its original shape and orientation. The corresponding SU(2) rotations, 
with the same axes and angles, are applied to point p in Fig. [44] taking us along a path 
p —> q — s- p' — >• p" . Point p" is not shown in Fig. 031 but it is a point of I\ that projects 
onto the same tetrahedron as point p, since the projected path in A±j is closed. Points 
p and p" differ by the phases of the four spinors, that is, by transformations generated 
by (Ji, I2, 13, 1 a)- Thus, there is a final segment p" —> p needed to close the path in 
, which runs along the T 4 fiber over the initial configuration in A4J . 

5.6. Closing the loop in $4^ 

There are several ways to compute the final four phases, but we will discuss just one. 
We start with vector Ji. The action of the rotations ([5"3"f on this vector can be written 

R(-U, 20 4 )i?(ji 2 , 20i 2 )Ji - Ji, (54) 

where we omit the subscripts on the i?'s because it is understood that only vector Ji 
is being acted upon, and where we omit the middle rotation in (I53[) since it does not 
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act on Jj. The product of the two rotations in (1541) is not the identity, but it is a 
rotation about axis ji since it leaves Ji invariant. 

To find the angle of this rotation, we use the Rodrigues-Hamilton formula 
(Whittaker 1960) for the product of two rotations in axis-angle form. Let a^, i = 1, 2, 3 
be three unit vectors, which we can plot on the unit sphere as in Fig. [48j We join 
the three points on the unit sphere by arcs of great circles. On following the path 
1 ^ 2 ^ 3 1 we regard the region to our right as the interior of the spherical 
triangle formed by the arcs. This gives meaning to the interior angles of the triangle, 
labeled (pi, i — 1,2,3 in the figure. On going from point i to point i + 1 we can go 
either the long way or short way around the great circle; the interior of the triangle 
and the definitions of the interior angles depend on which way we go, but the formula 
is valid in any case. If we follow the arcs the short way around, we obtain a spherical 
triangle such as that shown in Fig. UHJ Then the formula of Rodrigues and Hamilton 
is 

i?(a 3 , 2^ 3 )R{sl 2 , 20 2 )i?( ai , 20i) = /. (55) 

The proof is obtained by using ([48)) to write each rotation as a product of reflections, 
that is, i?(ai,2</>i) = Q(31)Q(12) and cyclic permutations, where for example Q(12) 
means reflection in the plane defined by axes ai and a 2 . 




Figure 48. Axes and angles 4>i, i = 1,2,3 for the product of three rotations. 

To apply (|55p to the rotations in (l54l) we imagine the vertex of the original 
tetrahedron (part (a) of Fig. |43|) where edges 1, 4 and 12 meet as at the center of the 
sphere in Fig. 051 and we identify axes (ai,a 2 ,a 3 ) with (ji 2 , — j4,ji). Then the 0's of 
Fig. 0FJ become the interior dihedral angles of the tetrahedron, and we have 

R(ji,24>M-U, 20 4 )i?(ji2, 20i 2 ) = I. (56) 

Thus, the product of the two rotations in (|5~4l is R(— ji, 20i). 

The third rotation found in this manner can be regarded as a holonomy. As the 
first two rotations in ([54)) are applied to Ji, that vector traces out a closed curve on 
S 2 which is composed of the arcs of two small circles (the curve is a "small lune"). 
The sphere in question can be regarded as Si, the symplectic manifold denoted £ in 
(|2~4l (the 1 subscript indicates that we are dealing with the first of the four copies of 
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S in £47). The two rotations themselves in (|54"f can be regarded as the lift of that 
curve on S 2 into 50(3). The lift is an open curve starting at the identity in 50(3) 
and ending at the product of the two rotations shown in (|54p . To close this curve in 
50(3), we apply the third rotation shown in (1561) . 

This closed curve in 50(3) may then be lifted (a second time) into SU(2), by 
replacing each 50(3) rotation by an SU (2) rotation with the same axis and angle, that 
is, the lift is specified by the product 2<fri)u(— J4, 204)it(jia, 20i 2 ). This product 
is either +1 or —1, depending on the homotopy class of the closed loop in 50(3). 
To find this class, we continuously deform the tetrahedron (see part (a) of Fig. |4"5)) , 
bringing dihedral angle 0i 2 to zero, so that the tetrahedron becomes flat. At the end 
of this deformation, 0i 2 = and one of tfii and </> 4 is and the other is n. Thus, the 
closed loop in 50(3) becomes an element of the noncontractible homotopy class of 
the homotopy group Z 2 of 50(3), so the product of the three SU(2) matrices is — 1. 
Therefore to close the loop in SU(2) we apply a final rotation by an angle — 2ir, to 
obtain 

u(h, -27r)u(h, 2<fc)u(-j 4 , 2<£ 4 )u(ji2, 2<f> 12 ) = +1. (57) 

The final rotation by angle —2tt could have been chosen about any axis if all we 
wanted to do was to close the loop in SU(2), but we choose axis ji for the following 
reason. The first two spin rotations in (|57p . when applied to a spinor in C 2 over 
vector J 1; produce another spinor that projects onto Ji again, that is, it differs from 
the initial spinor only by an overall phase. This overall phase is the U(l) holonomy 
mentioned above. The final step in closing the loop in Fig. 03] is to follow the / r -flows 
to adjust the phases of the four spinors, in particular, we must follow the /i-flow to 
adjust the phase of the first spinor. But the observable 3\ = I 2 is a function of I\, 
so we can follow its flow just as well. But J 2 generates a rotation about the direction 
Ji, that is, it is equivalent to multiplying the first spinor by an SU(2) transformation 
with axis ji and some angle. The angle required is seen in (f57|) : it is 2<pi — 2tt. 

We remark that the final rotation in (1571) could have been any angle — 2ir + Arm. 
Adding a multiple of 4-7T to this angle is equivalent to going around some closed contour 
on the A- manifold (in fact, within which, when the manifolds arc quantized and 
the Maslov phase is taken into account, changes the relative phase of the two branches 
by a multiple of 2tt. 

Next we turn to vector J3 and the phase needed to bring the third spinor back 
its original value after the rotations ((53)) . The action of these rotations on J3 is given 
by 

R(-U, 20 4 )i?(j 23 , 20 23 )J 3 = J 3 , (58) 

where we omit the first one since it does not act on J3. The product of the two 
rotations in (|58p is a third rotation about axis j 3 , which can be computed with the 
help of (|55[) . This time we identify the 3-4-23 vertex of part (c) of Fig. |45j seen more 
clearly in Fig. H7J with the origin of the sphere in Fig. [48] We use the version of the 
tetrahedron after the first rotation (part (c)) since the middle rotation in (|5"8")) involves 
the rotated axis j 23 . Thus we find 

i?(-j 23 , 2<Mi?(j 4 , 20 4 )i?(-j 3 , 20 3 ) = I. (59) 
Now using R(— a, a) = i?(a, a) -1 and taking the inverse of (f5"9"|). we obtain 

i?(j 3 , 2&)iZ(-j 4 , 204)^38.2^23) - /, (60) 
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which specifies a closed loop in 50(3). We find the homotopy class of this loop by 
taking 23 — > 0, which makes one of 04 and 03 zero and the other n, so the loop in 
5*0(3) belongs to the noncontractible class. Thus, the lift into 5(7(2) is not closed, 
but we can close it by appending a final spin rotation about axis J3 by angle — 2n. 
Thus the closed loop in SU (2) is specified by 

U(j 3 , -27l>(j 3 , 20 3 )u(-j 4 , 20 4 )u(j 23) 20 23 ) = +1, (61) 

which when applied to the third spinor at point p in Fig. |U returns it to its initial 
value. The final U (1) holonomy of spinor 3, after the application of the three rotations 
(PI), is 203 - 2tt. 

As for J 2 , its cycle on the sphere is specified by 

R(-H, 2fa)R(3' 23 , 20 23 )i?(ji2, 20i 2 )J 2 = J 2 (62) 

The product of the three rotations shown must be a rotation with axis j 2 . To find the 
angle, we first use to obtain 

i?(-j 4 , 204)^023,2023) = R(-j 3 , 203). (63) 

Substituting this into (|6"12|) we obtain a product that we can evaluate with the help of 
(|55p and with reference to the 2-3-12 vertex of the original tetrahedron (part (a) or 
(f) of Fig. US]). The result is 

i?(j 2 , 20 2 )i?(~j 3 , 20 3 )i?(ji 2 , 20 12 ) = I. (64) 

Thus the closed loop in 50(3) associated with the loop traced by J 2 on the sphere is 
specified by 

i?(j 2 , 2<h)R(-k, 204)i?(j 23 , 20 23 )i?(ji 2 , 20 12 ) = I. (65) 

To find the homotopy class of this loop we deform the tetrahedron into a planar shape 
as before, and find that two of the four angles (0 2 , 04, 23 , 0i 2 ) are and two are n. 
The loop in 50(3) thus becomes the product of two rotations with angles 2ir, which 
belongs to the contractible homotopy class. Therefore the lift into SU(2) is closed, 

ufa, 20 2 )it(-j 4 , 20 4 )u(j 23 , 20 23 )w(ji 2 , 20i 2 ) = 1, (66) 

which when applied to the second spinor at point p in Fig. [44] returns it to its initial 
value. The final U (1) holonomy of spinor 2, after the application of the three rotations 
US3J, is 20 2 . 

Finally, we treat vector J4 and the corresponding U(l) holonomy. The effect of 
(|53p on J4 is simply 

i?(-j 4 , 20 4 )J 4 = J 4 , (67) 

since the first two rotations do not act on J 4 . This specifies an open curve in 50(3) 
that can be closed (trivially) by multiplying by a rotation about axis J4 , 

R(h> 20 4 )-R(-j 4 , 20 4 ) = I. (68) 

The closed loop in 50(3) belongs to the contractible homotopy class, so its lift into 
SU(2) is closed and is specified by 

u(j4,20 4 M-j4,20 4 ) = 1. (69) 

When applied to the fourth spinor at point p in Fig. [U this sequence of spin rotations 
returns it to its initial value. The final U(l) holonomy of spinor 4, after the application 
of the three rotations (|53|) . is 204. 
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To summarize, we have succeeded in constructing the closed loop p-q-p'-p 
illustrated in Fig. as the product of a sequence of seven spin rotations, each one 
generated by the Hamiltonian flow of one of the observables in the A- or _B-list (in 
the A-list while we move on the A-manifold, and in the _B-list while we move on the 
J3-manifold). Most of these can be regarded as being generated by the square of some 
angular momentum vector; for example, the first and second spin rotations, specified 
by the axes and angles of the right-most rotations in (f53|) . are generated by Jf 2 and 
J23, respectively, while the last four rotations can be regarded as being generated by 
r = 1, . . . , 4. The third rotation in (|53")l , about axis j 4 , is generated by J4 • J, that 
is, u(— J4, 2^4) = exp(i^>4j4 • a) is applied to all four spinors. 

5.7. The Ponzano-Regge phase 

The actions associated with these spin rotations are easily computed, using the 
complex 1-form (|4"7} and the methods of Sec. 1.9.2. To summarize the results, let 
J p be a partial or total sum of the four angular momentum vectors (some range of 
r = 1, ... ,4 is summed over), with magnitude |J p | — J p . Then the action along the 
path generated by n • J p with elapsed angle 9 is simply (n • 3 P )9. In particular, the 
third rotation in ((55)) , the overall rotation of the tetrahedron taking us from part (e) to 
part (f ) in Fig. 23J does not contribute to the action since in this case J p = Jtot which 
vanishes on the A- and B-manifolds. See also (1.74). As for the rotations generated by 
some J^, in this case n = j p = J p /J p , so the action is simply J p 6. Thus the first and 
second spin rotations specified by (l53l) contribute Ji2(2</>i2) + ^23(2^23) to the total 
action. As for the last four rotations, in this case J p is one of the J r , r = 1, ... ,4, 
and the four angles are summarized in Sec. 15.61 These rotations therefore contribute 
Ji(20i — 2tt) + ^(2^2) + ^3(2^3 — 2tt) + J4(204) to the total action. Altogether, we 
have 

S=Ipdx = 2^2j r (j) r -2Tr(J 1 + J 3 ) 7 (70) 

•* r=l 

where index r = 5 means r — 12 and r — 6 means r — 23. This can be written 

S= -2* + 2tt(J 2 + Ji + J12 + J23), (71) 
where on quantized manifolds the final term is an integer multiple of 2-7r, and where 

6 

* = J2 M* ~ &•)■ ( 72 ) 

r=l 

The angle it — 4> r is the exterior dihedral angle, so VP is the phase of Ponzano and 
Regge. 

6. The amplitude determinant and reduced phase space 

Amplitude determinants are notorious for the trouble they cause in semiclassical 
approximations, for example, Gutzwiller's amplitude determinant (Gutzwiller 1967, 
1969, 1970, 1971) has that reputation and in several studies of asymptotic 
approximations to spin networks the authors have resorted to numerical calculations 
for the amplitude determinant. In fact, amplitude determinants can be expressed in 
terms of Poisson brackets, which aids considerably in their evaluation. For example, 
Wigner's (1959) amplitude for the 6j-symbol is a single Poisson bracket, while the 
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amplitude determinant for the 9 j-symbol is a 2 x 2 matrix of Poisson brackets (Haggard 
and Littlejohn 2010), the derivation of which was the easiest part in the asymptotic 
formula for the 9j'-symbol. Similarly it is easy to obtain tractable expressions for 
the amplitude determinant for the 15j-symbol and other more complicated cases of 
interest in quantum gravity. 

In Sec. 1. 11 we presented a coordinate-based discussion of amplitude determinants 
in the 3j-symbol. For a more geometrical treatment of some of the issues discussed 
there we refer to the literature on the "quantization commutes with reduction" 
theorems of Guillemin and Sternberg (1982). Here we will simply review the results 
of I and discuss their geometrical content. 

The discussion involves symplectic reduction, which in the case of the 4j-model of 
the 6j-symbol leads to the reduced phase space of the 6j-symbol, a 2-sphere denoted 
by r in (|25[) . This space is essential for understanding the semiclassical mechanics 
of the 6j-symbol, for example, it is the phase space that underlies the 1-dimensional 
WKB methods used by Schulten and Gordon (1975a, b), and it played an important 
role in the derivation of the uniform asymptotic approximation of Littlejohn and Yu 
(2009), as well as in the semiclassical studies of the volume operator in quantum 
gravity by Bianchi and Haggard (2011). 

6.1. Densities and amplitude determinants 

In this section we adopt a general notation, as in Sec. 14.31 so that our results can 
be applied to the 4j-model of the 6j-symbol, the 12j-model of the 6j-symbol, the 
9j-symbol, and many other examples. The phase space is P = M. 2N with symplectic 
form lo = dpAdx. The theory of symplectic reduction usually begins with a symplectic 
group action on phase space, but a more natural starting point in our applications 
is the components of the momentum map of the group, which is a list of classical 
observables. Actually, there are two lists. 

We assume that there exist on P an A- and a B-list of observables, (Ai, . . . , An) 
and (B\, . . . , -B/v), each of which forms a Lie algebra under the Poisson bracket. We 
denote contour values by a, and bi, so that level sets are specified by Ai = a.; and 
Bi = bi. We denote these level sets by L a and Lb (these are the A- and B- manifolds) . 
We assume that the A^s and B^s are functionally independent at most places in phase 
space, so that the generic dimensionality of L a and Lb is N, although it may be less 
in exceptional cases. The A- and _B-lists of observables correspond to two groups with 
symplectic actions on phase space, which we denote by G a and Gb, with Lie algebras 
Q a and Qb- We assume the groups are connected and have dimensionality N, as in 
Sec. 14.31 We denote the structure constants of one or the other group by c k ij, as in 
Sec. S3 

We are exclusively interested in contour values a and b that are fixed points of 
the coadjoint actions of the respective groups, since this implies that L a and Lb, if 
./V-dimensional, are Lagrangian. Only Lagrangian manifolds can support semiclassical 
wave functions. The restriction on the contour values means, however, that we do not 
have a Lagrangian foliation of P. 

We will be interested in the WKB or semiclassical x-space wave function 
associated with the A- or B-manifold. We will work with the A-manifold, since the 
B-case is similar. For modern treatments of WKB theory, see for example Martinez 
(2002) or Mishchenko et al (1990). The variables x are half of the coordinates (x,p) 
on P; we denote "x-space" by Q = M , which abstractly is best seen as the quotient 
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space when P is divided by the foliation into vertical Lagrangian planes, x = const. 
We write this wave function in the form 

1> a (x) = (x\a) = KY,Mx)\ 1/2 exp{t[S{x) ~ (73) 

br 

where if is a normalization; where the sum is over branches of the projection 
7r : L a — > Q from the Lagrangian manifold to Q, assumed to be locally invertible; 
where the branch index is suppressed in the sum; where S is the action computed as 
in Sec. 15.3) and where n is the Maslov index. 

We associate the function f2(x) with an ./V-form on Q by 

f2 = Q(x) dx 1 A ... A dx N . (74) 

It follows from WKB theory that f2(x) satisfies a set of amplitude-transport equations 
on Q, one for each observable Ai. These are conveniently expressed in terms of an 
iV-form a on L a , defined by a = 7r*f2. Then the amplitude transport equations on 
Q are equivalent to Cx i o~ = on L a , where C is the Lie derivative and where Xi, 
i = 1, . . . , N are the Hamiltonian vector fields on L a associated with the A+, that is, 
X.- L = oj~ 1 dAi. This means that a is invariant under the infinitesimal action of the 
group, hence under any finite action (recall that L a is an orbit of the group action). 
These are left actions. 

We express the solution a as follows. The vector fields Xi may be chosen as a 
(generally non-coordinate) basis in each tangent space at each point of L a . Let A 1 , 
i = 1, . . . , N be the dual 1-forms on L a , that is, X t (Xj) = <5j. Then, in the case of 
compact groups, the solution a is given by 

(7 = A 1 A ... A X N , (75) 

to within a multiplicative constant. The reason is that a is a version of the right- 
invariant Haar measure on the group, and for compact groups the left and right Haar 
measures are equal. Thus, a is invariant under the (left) action of the group, and 
C Xt a = 0. 

The action of G a on L a provides two ways in which geometric structures on 
the group can be transferred to L a . First, the group action implies a linear map 
■ 8a — > T x L a for each x G L a , which, under our assumptions, is invertible. This map 
can be used to push forward a standard X-form on Q a to L a . Let the basis in Q a be 
= 1, . . . , N}, corresponding to the observables Ai and vector fields Xi, and let 
{a l ,i = 1, . . . ,N} be the dual basis in g*, that is, a l (^j) = Sj. Then if the iV-form 
a 1 A ... A a" is pushed forward to L a in this manner, we obtain a, defined by (|75|) . 

A second way involves picking a point xo € L a to serve as an "origin" in L a , and 
then identifying points x of L a by the group element g such that x = gxo. This creates 
a diffeomorphism between a neighborhood of the identity in G a and a neighborhood of 
xq in L a , which can be used to push forward differential geometric structures from G a 
to L a . It then turns out that a given by (|75[) is the push forward of the right- invariant 
N-iorm on G a associated with a 1 A ... A a . Thus, a is naturally invariant only under 
a right action of the group; the only reason it is also invariant under the left action is 
that for the groups we consider, the left and right Haar measures are identical. 

To see this from another standpoint, a short calculation shows that Cx^ = 
c : 'ik\ k , from which follows Hxi°~ — c'ij To derive this we recall that although 
the Poisson brackets {A^A,} vanish on L a , the Lie brackets of the corresponding 
Hamiltonian vector fields do not, instead we have [X U XA = ~X k c k , Lj . But if the left 
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and right Haar measures are equal, then the adjoint representation Ad g is volume- 
preserving, so the structure constants are traceless, and C,Xi a = 0. In our examples 
we deal only with compact groups, so this condition is met. 

That the solution is unique to within a multiplicative constant can be seen by 
supposing that a' is another TV-form on L a such that Cx^ 1 — 0. Since all N- forms 
are proportional, we must have a' = fa where / is a function on L a . This implies 
x i(f) = 0, or / = const on L a . 

Now given that (|75p is the solution we want, we write a = ir*Q in the form, 

Q(x) dx 1 A ... A dx N = A 1 A ... A X N , (76) 

where we write simply dx 1 for Tr*(dx l ), and we evaluate both sides on the set of vectors 
(Xi, . . . ,X N ). This gives 

Q{x)det{x i ,A j } = 1, (77) 

where we use dx l (Xj) = {x % ,Aj}. This reproduces Eq. (1.92). 

Finally, as shown in I, the normalization integral, evaluated in the stationary 
phase approximation, implies K = 1/ y/Va, where V a is the volume of L a with respect 
to a. 

6.2. The amplitude of (b\a) 

A coordinate-based derivation of the amplitude of the semiclassical matrix element 
(b | a) was presented in I for the case of the 3j-symbol. Here we discuss the results from 
a geometrical point of view, using the general notation of Sec. 14.31 and the previous 
section. The main issue is that the stationary phase set L a l~l Lb in the evaluation of 
the integral (b\a) — j dx(b\x){x\a) is not a set of isolated points, as expected on the 
basis of a naive dimensionality count, but rather a set of manifolds of dimensionality 
> 1. These nontrivial intersection manifolds are due to the existence of a common 
"intersection group" of G a and Gb, which may be defined as follows. 

The basic idea is that the A- and £>-lists of observables may have some observables 
in common, which generate the intersection group. This is obviously the case in (|27j) . 
for example. But the specific observables that occur in the A- and £>-lists depend on 
the bases chosen in the Lie algebras g a and g&, and if a basis is changed then the A- or 
£>-observables are replaced by linear combinations of themselves. So we need a precise 
definition of the observables "in common." 

The actions of groups G a and Gb on P provide Lie algebra anti-homomorphisms 
between g a and g& and the Lie algebra of (globally) Hamiltonian vector fields on P. 
By our assumptions these anti-homomorphisms have full rank, so the images of these 
two anti-homomorphisms are two iV-dimensional Lie algebras of Hamiltonian vector 
fields. These two Lie algebras have an intersection which itself is a Lie algebra. Let 
the dimension of the intersection be p, and let p + q = N. The intersection Lie algebra 
is generated by a set of Hamiltonian functions, call them (C%, . . . , G p ). These can be 
regarded as the functions common to the original A- and B-lists, that is, by a change 
of basis in Q a and Qb we can bring the A- and _B-lists into the form A = (C, D) and 
B = (C, E), where D = {D\, . . . , D q ) and E — (Ei, . . . , E q ) are sets of observables 
that the A- and B-lists do not have in common. The common observables C generate 
the action of an "intersection group" Go, that is, they are the components of the 
momentum map of the action of Go on P. The group Go is not uniquely determined 
by the G observables, only its action on P is. The same applies to G a and Gb, which 
are generated by the A- and B-lists of observables. But in practice there are convenient 
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ways of choosing all these groups so that the kernels of their actions are discrete. This 
means that dimG a = dimGb = N and dim Go = p. We may assume moreover that 
Go is connected (as we have already for G a and Gf,), since we are free to take the 
connected identity component of any of these groups. 

The level set of the momentum map of the intersection group Go plays an 
important role in what follows. We denote this level set by Lq; its equation is d = a, 
i = 1, . . . ,p, for some contour values ct. We choose the Ci such that Lo has its generic 
dimension 2N — p = N + q > N . 

The level set Lo is foliated into A- and S-manifolds, where the A-manifolds are 
parameterized by the values di of the observables Di, and the B-manifolds by the 
values ej of the observables Ei, i = 1, . . . , q. We will assume that the Poisson brackets 
{Ai, Aj} and {B i: Bj} vanish on L , so that the generic A- and _B-manifolds in L are 
Lagrangian. This is an extension of our earlier assumption, that these Poisson brackets 
vanish on a specific pair of A- and _B-manifolds. Thus, after removing exceptional 
manifolds of less than generic dimensionality, Lq is foliated into Lagrangian manifolds 
in two different ways. These Poisson bracket relations imply {Ci,Cj} — on Lq, so 
that c is a fixed point of the coadjoint action of Go- This does not mean that Lo 
is Lagrangian (the dimension is wrong), but it is coisotropic (Abraham and Marsden 
1978). It also means that the isotropy subgroup of the coadjoint action is the whole 
group, so the reduced phase space of level set Lq under the Go action is the space 

This draws attention to the orbits of Go, which generically have dimension p. 
Not only is L foliated into orbits of Go, so is each A- and B- manifold, since Go is 
a simultaneous subgroup of G a and G&. Thus, the intersections of L a and Lb, which 
are the stationary phase sets, are also foliated into orbits of Go- We will assume that 
L a n Lb is a union of a discrete set of orbits of Go- Based on a dimensionality count, 
this is the generic case. It holds for example in the 4j-model of the 6j-symbol, where 
L a n Lb consists of two orbits of Go, the sets I\ and 1% in Fig. HH It also holds in the 
12j-model used by Roberts (1999) and in our own work on the 9j-symbol (Haggard 
and Littlejohn 2010). 

The orbits of Go appear in the stationary phase evaluation of 

(b\a) = [ dx (b\x)(x\a) = . / dx~S~^ — . 

Xl/ J VIAI/ jVaVbJ ^ ^\&ct{x,A}dct{x,B}\ 

xexp{i[S A (x)-S B (x)- fJ ,Tr/2}}, (78) 

where we have inserted the normalized x-space wave functions for the states \a) and 
\b), where the sum is over all pairs of branches of both wave functions, and where fj, 
is the cumulative Maslov index. The stationary phase set is the projection of L a n Lb 
onto Q (cc-space); it is a union of the projections of the discrete set of orbits of Go 
that make up L a C\ Lb- 

The projection of each orbit is a subset of Q that is locally p-dimensional. As 
in I, we introduce a local change of coordinates x —¥ (y,z), where (y , . . . , y p ) are 
coordinates along the projected orbits and (z 1 , . . . , z 9 ) are coordinates transverse, 
with z — being the projected orbit itself. Then the integral becomes (suppressing 
normalization and branch sums) 

— / g-L exp{z[S A (y, z) - S B (y, z) - ^/2]}. (79) 

det{y,C}\J y/\det{z,D}det{z,E}\ 

When the exponent is expanded to second order in z, the leading term Sa(v, 0) — 
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Sb(u, 0) is independent of y, that is, it is constant along the orbits of Go- The second 
derivative matrix of Sa — Sb with respect to z does depend on y, but after doing 
the Gaussian integral and combining with the determinants in the denominator of the 
^-integration, the result is expressed purely in terms of the Poisson brackets {Ei, Dj}, 
which are independent of y. 

The fact that the z-integral is independent of y, that is, our location on the orbit 
of Go, is noteworthy. It means that the integral does not depend on the detailed nature 
of the z-coordinates, for example, it is invariant under a coordinate transformation of 
the form z' — z'(z, y) such that 2 = implies z' = 0. It suggests that we are dealing 
with a quotient operation in which we divide by the orbits of Go- 

Another remark is that there is nothing special about the x-representation in 
which the integral (|75|) is carried out. The ^-coordinates were introduced as half of 
the (x,p) coordinates on P — M. 2N , but any representation related to this one by a 
metaplectic transformation (Littlejohn 1986) would work as well. This amounts to 
foliating phase space, not by the vertical Lagrangian planes x = const, but rather by 
other Lagrangian planes related to this one by any linear, symplectic map. In this 
manner one can divide the orbits of Go into segments and do the integral over each 
segment in a representation in which the projection of the orbit onto the representation 
space has full rank. The segments into which the orbit of Go is divided can even be 
infinitesimal, effectively making the representation for the integral a function of where 
we are along those orbits. Similarly, one can avoid the caustics of the wave functions 
ip a and ipb- This is the old idea underlying the Maslov method (Maslov and Fedoriuk 
1981) in WKB theory. 

In any case, once the z-integral is done and is recognized to be independent of y, 
the remaining y-integral can be lifted to the orbit of Go in phase space whereupon it 
becomes just the integral of the Haar measure of Go, giving the volume of the orbit. 
This Haar measure is normalized as in the previous section, that is, we start with 
the observables Ci, i = l,...,p, we associate these with Hamiltonian vector fields 
Xi = Lu~ 1 dCi (a change of notation from above, where the X's were associated with 
the A's), we define form A 1 dual to the Xi on the orbits of Go, and then the Haar 
measure on the orbit is taken to be A 1 A ... A A p . The final result is 



where Vj is the volume of intersection manifold / (an orbit of Go), where Si means 
Sa ~ Sb evaluated on intersection manifold /, and where fi is a cumulative Maslov 
index (not necessarily the same as the previous ones). The branch index could 
otherwise be written as /, a label of the intersection manifold, and Vj is taken out 
of the sum because it docs not depend on which intersection manifold is taken. The 
volume Vj differs from the volume of Go because in general there is a discrete isotropy 
subgroup, as in (1441) (one is really computing the volume of a coset space). 

The result (|80|) contains a q x q matrix of Poisson brackets, {Ei,Dj}, whose 
geometrical content may be understood in terms of a variation of the discussion of 
densities in Sec. 16.11 First we recall that the A- and B-lists are decomposed according 
to A — (G, D), B = (G, E). Next, we fix the contour values Cj, i = 1, . . . ,p, so that 
we have a definite level set Lo of the momentum map of Go- Then we let "&-space" be 
ffi 9 with coordinates (ei, . . . , e q ) or the region of M. q that is the projection of Lq onto 
K 9 , where coordinates are interpreted as the contour values in Ei — ei. The matrix 




(80) 



where now the branches refer to the discrete set of orbits of Go that make up L a C\Lf,, 
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element (b\a) can be thought of as a wave function on 6-space for fixed values of the a's, 
that is, of the c's and d's. We write the amplitude of the semiclassical approximation 
to this wave function as 117(e)! 1 / 2 , where fi = £l(e)dei A ... A de q is the associated 
density (a g-form on &-space). This density is the projection of the natural density on 
the A-manifold, in the following sense. Let the A- or (G, D)-observables be associated 
with vector fields A,; = us~ l dd : i = 1, . . . ,p, and Yi = us~ 1 dDi, i = 1, . . . , q, with dual 
1-forms A\ i = l,...,p and i = 1, . . . , q. These induce a density ctq A fi 1 A . . . A \i q 
on L a , where <7o = A 1 A . . . A X p is the Haar measure on Go- Then fl satisfies 

(T A TT*n = (T A fjt 1 A . . . A fi q , (81) 
where n is the projection from Lq or L a C Lq onto 6-space. Now evaluating both sides 
on the set of vectors (X\, . . . , X p , Yi, . . . , Y q ), we obtain 

fl(e) det dE t (Yj) = 1. (82) 
This gives the amplitude shown in since dEi(Yj) — {Ei,Dj}. 

This discussion has treated the A- and B-manifolds asymmetrically, projecting 
from the A-manifold onto 6-space, but we could have projected the density on the 
B-manifold onto a-space with the same result. There are really four densities, two on 
the A- and B- manifolds, and two of the A- and B- manifolds. 

This discussion leads us to consider the reduced phase space T = L /G , which 
is parameterized by the contour values Cj, i = 1, . . . ,p. As is standard in symplectic 
reduction, the symplectic form on T is obtained by pulling back vectors from T to 
Lo C P and evaluating them on the symplectic form on P; this is meaningful because 
the answer does not depend on where on the orbit of Go they are pulled back to, nor 
on the component of the pulled-back vectors along the orbit. A consequence is that 
the projections of the A- and P-manifolds onto T, which are g-dimensional since the 
Go orbits are p-dimensional, are Lagrangian on T. Since we are assuming that L a n L 
is a discrete union of Go orbits, the projected manifolds on T intersect in a discrete 
set of isolated points. This would be the generic case on a symplectic manifold of 
dimension lq. Also, Go-invariant functions on Lq project onto functions on T, whose 
Poisson brackets on T are the same as the Poisson brackets of the original functions 
on P. Such functions in the present discussion include the observables Di and Ei, 
i = l,...,q, so the Poisson brackets {Ei,Dj} of ([80|) are naturally interpreted as 
living on T. 

6.3. The case of the 6j -symbol 

It is straightforward to apply the general theory of Sees, l6.l land l6.2l to the case of the 
6j-symbol. The phase space is C 8 = M 16 so N — 8. The common observables G are 
J r , r = 1, . . . , 4, and Jtot, so p — 7 and q = 1. The group Go is T 4 x 51/(2). The level 
set Lq of the momentum map of Go, for = J r , r = 1, ... ,4, and J to t = 0, is the 
subset of <&4j upon which the four angular momenta have specified lengths and their 
vector sum is zero. It is logical that this would be the subset of the classical phase 
space $4j that corresponds to the space Z of four-valent intertwiners, introduced in 
Sec. 13.41 on which I r = j r , r = 1, . . . , 4 and J to t = 0. 

The contour values of the G's, that is, of I r , r = 1, . . . , 4 and Jtot, must be chosen 
so that Lq has its maximum dimensionality, namely, 9. The condition is J r > 0, 
r = 1 , . . . , 4 and the polygon inequality, 

l 4 

max{ Ji, J 2 , J 3 , Ji) < - ^ Jr > ( 83 ) 

r=l 
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which as indicated must not be saturated. This is the condition that it is possible to 
make a noncollinear polygon in ]R 3 out of vectors of four given positive lengths. If this 
condition is satisfied, then Lq is foliated into A- and B-submanifolds whose generic 
dimensionality is 8. 

There is only one observable of the D- and E'-type; according to (|27| we should 
make the identifications D — J 2 2 ancl E — J23. However, since we used d9 in the 
volume form on the A- and i?-manifolds when computing the volume (|34[) , we should 
use instead D = | J12 1 and E — | J23], since these are conjugate to 9 (really 9 a and 9b, 
since there are two of them) . Then the Poisson bracket for the amplitude is 

{E, D\ _ { |J 2 + ,,|, |J, + ,,|) . hJh+g£!h±M . " m 

I J 12 1 |J 23 I 

where V is the signed volume of the tetrahedron, 6V = J\ ■ (J2 x J3), where we 
have used ff^B]) to evaluate the Poisson bracket, and where in the final step we have 
evaluated the Poisson bracket on A- and B-manifolds with contour values J12 and 
J23. We see the appearance of Wigner 's volume (Wigner 1959). The volume changes 
sign between the two intersection manifolds I\ and I2 , but it appears with an absolute 
value sign in (|80p so both stationary phase points in the 6_7-symbol have the same 
amplitude. 

The amplitude If^ 1 / 2 contains the factor 

1 



y/j^hs = - v/(2j 12 + 1)(2j 2 3 + 1), (85) 

which, when evaluated as shown on quantized manifolds, reproduces the square roots 
seen in (|19j) . Thus, based on the pieces of the formula we have determined so far, we 
can write 

e W 4 1 



n ]2 J12 



e i(Si-/W2) _j_ e i(S 2 -M27r/2) 



. . , . (86) 

33 34 323 J y / l2n\V\ 2 

where Si and S2 are given by (|45|) (they are the phases Sa — Sb, evaluated on 
the intersection manifolds I\ and I2 of Fig. |43|) . These phases contain the phase 
conventions for the states \A) and \B) shown in Fig. [37l which appear semiclassically 
as the origins for action integrals on the Lagrangian manifolds, denoted ao and 60 in 
Fig. 031 Since we have not considered these phase conventions yet, we cannot say what 
Si and S2 are, but the difference S = S2 — Si is given in terms of the Ponzano-Regge 
phase by (|7T|) . 

To return to the 9-dimensional space L , it projects onto what was called "6- 
space" in Sec. 16. 2[ which in the present case is the J23-axis, producing the interval 
[>^23,min: <^23,max], given by (H?2"j) . The inverse image of a point on this interval is a 
£?-manifold with a given value of J23. Similarly Lq projects onto "a-space", that is, 
the Ji2-axis, producing the interval (1311) . the inverse image of a point of which is 
an A-manifold with given J12 value. The space Lq projects onto the J12-J23 plane 
producing the classically allowed region of Fig. |4T] the inverse image of a point of this 
region is an intersection of an A- and a £>-manifold of given J12 and J23 values. If the 
point does not lie on the caustic curve, the intersection consists of two disconnected 
components Ii and I2, each an orbit of Go, as in Fig. [43] On the caustic curve these 
components merge into one. Finally, by dividing by Go, Lq projects onto the reduced 
phase space T, which has dimensionality 9—7 = 2. We now turn to this space. 
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6.4- The reduced phase space T 

Spaces of the type T seem to have appeared first in the work of Kapovich and Millson 
(1995, 1996). Those authors showed that the space of polygons of a given number of 
sides with fixed lengths in R 3 , modulo overall rotations, is a symplectic manifold. In 
fact, for quadrilaterals the space is precisely T. This space was recently subjected to 
direct geometric quantization by Charles (2008), who connected it with the 6j-symbol 
and used it for a new derivation of the Ponzano-Regge asymptotic formula. The space 
of five-sided polygons is the analog of T for the 9j-symbol; it was used by Haggard 
and Littlejohn (2010) in their study of the asymptotics of the 9j-symbol. 

To visualize T = Lq/Gq it helps to carry out the reduction in two steps. In 
the first step, we choose contour values J r > 0, r = 1, . . . , 4, for which the level set 
= J r in <E>4j is the product of 3-spheres (S 3 ) 4 , as shown in (|25|) . The I r ,r = 1, . . . , 4, 
generate the action of the group U(l) 4 , corresponding to the phases of the four spinors. 
Dividing this level set by U(l) 4 , we obtain the symplectic manifold £4^ shown in ([231) . 
which topologically is (5 2 ) 4 , and which consists of all sets of four vectors 3 r in R 3 , 
r = 1,...,4, with fixed lengths, |J r | = J r . The space £4^ is 8-dimensional. In the 
second step, we consider the submanifold in £4^ upon which J tot = 0, which is a level 
set of the momentum map of the action of 5*0(3) on £4^. This manifold consists of 
sets of four vectors 3 r in R 3 of fixed lengths J r such that ^ J r = 0. The vectors 
can be placed end-to-end to form a "closed link," that is, a four-sided polygon in R 3 . 
The set of closed links is denoted CL in (f25|) . We assume the polygon inequality (f83| 
is satisfied, so the space CL has dimensionality 8 — 3 = 5. Since this is the level set 
Jtot = 0, the isotropy subgroup of the S f O(3)-action is 5*0(3) itself, so the reduced 
phase space is CL/S'0(3), which is the space T in (l2"5j) . This has 5 — 3 = 2 dimensions. 

The phase space T is parameterized by the four fixed, positive values J r , 
r = 1, ... ,4 that satisfy the polygon inequality (|83)l . A point of this space specifies 
a quadrilateral in R 3 of the given lengths, modulo overall proper rotations. Once the 
quadrilateral has been determined, one can draw in the two remaining edges, of lengths 
J12 and J23, to fill in a tetrahedron. Thus, T can be thought of as the shape space for 
a set of tetrahedra, four of whose edges have fixed, positive lengths. The two lengths 
that are variable are on opposite sides of the tetrahedron. Here we define "shape" 
as a configuration modulo proper rotations, as in Littlejohn and Reinsch (1997); two 
shapes related by spatial inversion are generally distinct. The lengths | J12 1 and | J23 1 
are variable functions on T. In fact, any rotationally invariant quantity associated 
with the tetrahedron, such as the dihedral angles, the areas of the faces, the signed 
volume, etc, is also a function on L. 

It is easy to find coordinates on L. We may take one coordinate to be J12, which 
varies between the bounds (|31|) . (From this point on we drop the distinction between 
IJ12I and J12, and similarly for J23 and J13.) For a fixed value of J12, the allowed 
set of shapes is generated by executing the "butterfly" motion about the axis J 12, 
that is, rotating the triangle 1-2-12 about this axis, relative to the 3-4-12 triangle. 
We recall this motion is the Hamiltonian flow generated by Jf 2 or the magnitude 
I J12 1 (see (|33|)). Thus, coordinates can be taken to be (Ji2,0i2), where 0i2 is the 
interior dihedral angle about the 12-edge. For each value of J 12 in the interior of 
the range (j3"Tj) . a circle of shapes is generated as ^12 goes from to 27r; but at the 
endpoints there is only a single shape. For example, at the lower limit of (j3"Tj) . the 
case Ji2,min = | Ji — J2 1 is illustrated in part (a) of Fig. [49] In this case the rotation 
of vectors Ji and J2 about the axis J12 does not change the shape. Similarly, part 
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(b) illustrates the case J 12 = Ji2,min — 1^3 — Ji 1- In this case, the J^-action rotates 
the 1-2-12 triangle, but does not change the shape since the new configurations that 
result are related to the original ones by an overall SO(3) transformation. A similar 
analysis applies at the upper limit of f|3X|) . 




J 1 J 4 



Figure 49. When J12 is at its lower limit for given ( Ji, J2, J3, J4), the shape of 
the tetrahedron is not changed under the jJ 2 -action, so the lower limit corresponds 
to only one point of the reduced phase space. 

We see that the set of shapes generated in this manner for J12 in the interior of 
the interval pip is a cylinder, but the two endpoints are single points that pinch the 
cylinder at the two ends, creating a topological sphere. Topologically, Y = S 2 . 

The symplectic form on T may be obtained by projecting lj = dp A dx on <&4j 
as described in Sec. 16.21 , but it is easier just to notice that </>i2 is the parameter of 
evolution along the flow generated by J12, so (^12, J12) form a canonically conjugate 
(q,p) pair on the sphere. The same obviously applies to J23 and J31 and their conjugate 
(dihedral) angles, so we have 

dJ 12 A d(f> 12 = dJ 23 A d(p23 = dJ 31 A #31, (87) 

indicating three sets of canonical coordinates on T, related by canonical 
transformations. These are examples of the action-angle variables discovered by 
Kapovich and Millson (1995, 1996), which in all cases are closely related to the 
recoupling schemes used in angular momentum theory. The length J31 and the 
associated dihedral angle 03i do not appear in the coupling scheme we described in 
Sec. 13.41 or in the tetrahedra we have discussed so far, but would appear in a different 
tetrahedron in which the four vectors are placed end-to-end in a different order. 

So far we have described the construction of Y as a purely classical problem, but 
if the J r are quantized, J r = j r + 1/2, r = 1, ... 4, then one can speak of a quantized 
level set Lq and quotient space T. 

The area of the sphere T with respect to the form dJ\ 2 A d<p\2 is obviously just 
27r(Ji2 !max — Ji2,min)- If T is quantized, then by (|38[) and (|T6|) the area is where 
D = dim Z. When quantized, Y contains one Planck cell of area 2tt for every state in 
the Hilbert space Z. Obviously we obtain the same area if we use either of the other 
symplectic forms in (|87[). 

6.5. Quantized curves in Y 

States in Z, such as the A- and B-states given by (|2"Tj) . are represented semiclassically 
by Lagrangian manifolds in T, which are quantized curves in that space. For example, 
the A-states are represented by quantized levels sets of J\ 2 . 

To plot these we map Y into a unit 2-sphere in R 3 with standard coordinates 
(x, y, z) by associating (</>i2, J12) with a standard set of spherical angles (9, </>), where 
<p = 0i2 and 

Jl2 = Jl2,mm + y(l + COS A), (88) 
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where D = Ji2. m ax — ^i2,min (generally), or D — dimZ (when T is quantized). Here 
x = sin 9 cos <fi, y = sin 9 sin cf>, z = cos 9 are understood. Then the symplectic form on 
the sphere is (—l/2)Dsm9 d9 A dip. These coordinates make the Jyi orbits look nice 
(they are small circles), but not the orbits of the other variables, so they should not 
be thought of as having any privileged role. We use them mainly for plotting figures. 
The x-z plane in these coordinates does have an invariant meaning, however; it is the 
plane upon which the tetrahedra are planar (the volume vanishes). Moreover, time 
reversal is a reflection in this plane (y —> —y). 

The quantized orbits of J12, on a space T with five Planck cells of area, are 
illustrated in part (a) of Fig. [50] The orbits are just small circles, as noted. The orbits 
are numbered n = 0, . . . , 4, in order of increasing J12, and enclose area (n+i)27r. The 
minimum and maximum values of J12 are at the south and north poles, respectively, 
because of the choice (|88l) of coordinates, but the minimum and maximum quantized 
values are separated from the values at the poles by i, as indicated by (159"1) . since the 
polar caps defined by the last quantized orbit before the poles must have area \. 

The quantized orbits of J23 are illustrated in parts (b) and (c) of that figure. 
The function J23 has a minimum on the x-z plane with x > 0, and increases 
monotonically toward a maximum on the same plane with x < 0. The orbits are 
time reversal symmetric, but have no symmetry under inversion through the origin 
in the coordinates used (there is no reason why they should). Also shown are some 
quantized orbits of J13, in part (d) of Fig. [50] 

Another interesting observable on the reduced phase space is the signed volume 
V = (l/6)Ji • (J2 x J 3 ). This has been suggested by Chakrabarti (1964) and by Levy- 
Leblond and Levy-Nahas (1965) as a "democratic" alternative to the usual choices Jf 2 , 
J| 3 , J^ 3 for the intermediate variable in the quantum problem of the recoupling of three 
angular momenta. More recently, Carbone et al (2002) have reconsidered the use of 
this observable, and have derived both recursion relations and asymptotic formulas for 
the matrix elements connecting the Ji2-basis and the F-basis. The level sets (orbits) 
of V are plotted in part (e) of Fig. [50] The orbits plotted are not quantized, but are 
evenly spaced in V about V = 0, and labelled n — —2, . . . , +2. The V — contour is 
the x-z plane, which is also a quantized orbit when dim Z is odd. The orbit V = v is 
mapped by time reversal into the orbit V = —v. 

The Bohr-Sommerfeld condition for the quantized orbits of V can be expressed 
in terms of complete elliptic integrals of third kind. We will report on this and other 
results on the spectrum of V and on the wave functions (ji2|V") in future publications. 
Preliminary results are reported in Bianchi and Haggard (2011). 

6.6. The 6j -symbol on the reduced phase space 

As noted below (fT9|) the 6j-symbol is proportional to the matrix element (J23IJ12), 
so on r it is represented semiclassically by the quantized curves, J23 = j'23 + \ an d 
•hi = 3\i + \- This is illustrated in part (f) of Fig. [50] for the case (j'23 IJ12) = (0|4) , 
where the numbers refer to the labeling of the quantized curves in parts (a) and (b) 
of that figure. These curves are projections of the A- and _B-manifolds in $4^, and 
their intersections, the points labelled P and Q in the figure, are the projections of 
intersections manifolds I\ and I2 as illustrated in Fig. [43] respectively. This follows 
because we defined 1\ as the intersection manifold upon which V < 0, and V < at 
point P since it lies in the region y < 0. Also, the path (with direction indicated) 
surrounding the shaded area in part (f ) of Fig. [50] is the projection of the loop (with 
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Figure 50. Orbits (or level sets) on the reduced phase space. In (a), orbits of 
J12; in (b) and (c), of J23; in (d), of J13; in (e), of the volume V. Part (f) shows 
the orbits relevant to the 6j'-symbol. 

direction indicated) in $4., illustrated in Fig. [44] We assign arrows to the path in T 
by following the flows generated by Jf 2 and J| 3 in the direction of increasing </>i2 and 
<p23- Thus, the path The path P — > Q — > P in part (f) of Fig. l50l corresponds to the 
sequence of rotations carried out on the tetrahedron in Fig. 05l On reaching point 
Q, the path turns down, because the rate of change of J12 along the J23-flow is the 
Poisson bracket { J12, ^23}, which is negative when V > 0, as follows from (|84[) . One 
can also see this geometrically in Fig. 25] parts (b) and (c) of that figure correspond 
to point Q, and it is clear that J12 at first decreases, then increases again, on carrying 
out the rotation about J23, that is, on passing from part (c) to part (d) of Fig. 03 
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The symplectic area of the shaded region in part (f ) of Fig. [SU] is the same as 
the Ponzano-Regge phase, to within a constant that depends on j r , r = 1, . . . , 4. On 
quantized manifolds it will not matter what region on S 2 is taken as the interior of 
the closed loop on the sphere (the choices differ by multiples of the total area of the 
sphere), but to be precise one must worry about this. The Ponzano-Regge phase is 
also —2 times the i^-type generating function (Goldstein 1980, Miller 1974). Let us 
write (q,p) = (0i2, J12), (Q,P) = (023, J23) and F 4 (p,P) = -(1/2)5", where S is given 
by fTIjl. Then we have 

'--w-Sat**- m 

r— 1 

where we use index 5 instead of 12. Here the angles <fi r must be interpreted as functions 
of the J's. 

A more geometrical way of stating the same thing is to express the canonical 
transformation (0i2, J12) —> (023, ^23) as a map between two spheres with symplectic 
forms (|87[) ■ Then the graph of the canonical transformation is a Lagrangian manifold 
in S" 2 x S 2 with symplectic form <iJ23 A d023 — dJi2 A d0i2- This Lagrangian manifold 
is itself a sphere, whose projection onto J12-J23 space is precisely the interior plus 
boundary of the oval curve in Fig. 1411 The Lagrangian manifold has a double- valued 
projection onto the interior of the oval region, and a single-valued projection at the 
boundary (the caustics). The two disks fit together to form a sphere. 



6. 7. The Ponzano-Regge Formula 
If we factor out the phase 

' Sl+ S 2 \ { fil + H% \ 7T 



2 ™ 



from the quantity in the square bracket in (|86[) . we obtain a quantity proportional to 

(91) 



cos 



- (Ml ~ t*2)j 



However, in view of (|4"5l) . (|4*51) . (jTTj) and (fT2")l . this can be written as a sign times 



7T 



cos(* + A|U-J, (92) 

where the sign is (— Ipa+M+Jia+jaa anc j where A/x = H2 — Mi- Thus to obtain the 
argument of the cosine we must compute the relative Maslov index between the two 
branches, just as we have already computed the relative action S — S2 — S%. 

This is the easiest of the Maslov index computations in deriving the Ponzano- 
Regged formula. We will omit details and simply remark that the calculation can be 
carried out entirely on the reduced phase space T, where it is not difficult to show 
that Afi = 1. 

There remains the overall phase, which, in view of the reality of 6j-symbol, must 
be a sign ±1. This phase cannot be computed without taking into account the phase 
conventions for the two states \A) and \B), as well as the absolute Maslov indices on 
the intersection manifolds I\ and I-i- The analogous phase for the 9j-symbol was the 
most difficult part of the derivation of the results presented in Haggard and Littlejohn 
(2010). Since (in the case of the 6j-symbol) the answer is known, we can see that this 
overall phase (combined with the e l7T / 4 appearing in (|86]l ) must be 1. We had no such 
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luxury in the case of the 9j-symbol, nor did we find it possible to guess the answer. 
In any case, the final result is the Ponzano-Regge formula, 



For reasons of space we will omit the details in the calculation of these final 
phases, also because we have made the main points we wanted to make about the 
6j-symbol. These are the importance of the reduced phase space T, the geometry of 
symplectic reduction connecting it with higher dimensional spaces, and the manner in 
which elements of the semiclassical calculation (phases, amplitude determinants, etc) 
can be mapped from one space to another. 

7. Conclusions 

A glance at the calculation of Roberts (1999) shows that it is much easier and more 
elegant to compute the relative phase S = S2 — S\ of the 6j-symbol in a symmetrical 
or 12j-model than in the 4j-model presented in this paper. The sum over edges times 
dihedral angles appears almost immediately. On the other hand, we have given a 
much easier way of computing the amplitude determinant, reducing it to a single 
Poisson bracket, whereas Roberts had to evaluate a large determinant. Obviously 
what is needed is a way of connecting the various models at a semiclassical level, so 
that actions, amplitude determinants, Maslov indices, etc, can be mapped from one 
model to another and computed wherever most convenient. We have investigated this 
question and will report on our results in the future. For now we merely make a few 
comments. 

Already in this paper we have mapped amplitude determinants between various 
spaces, such as the Schwinger phase space $4^ and the reduced phase space T, which 
are connected by symplecic reduction. One might guess, therefore, that the 12j- or 
symmetrical model of the 6j-symbol and the 4j-model are related by some kind of 
symplectic reduction. This is not the case, however, at least not in the way that <&4j 
and T are related. One piece of picture, however, is the following. 

It is well known that a unitary map on a Hilbert space U : % — > % corresponds 
semiclassically to a symplectic map or canonical transformation T : $ — > <!>, where 
$ is the symplectic manifold corresponding to H. Also, the symplectic map T is 
conveniently viewed via its graph in $ x $, regarded as a symplectic manifold in its 
own right with symplectic form u>i — u>2, where the subscripts refer to the first and 
second factors of $ x $ (Abraham and Marsden 1978). With this understanding, the 
graph of T is a Lagrangian manifold, one which supports semiclassically the operator 
U in the same way as a Lagrangian manifold in $ supports a vector in %. That is, 
U, which begins as a map : % — > %, is reinterpreted as a function : H <8> H* — > C, 
that is, a "wave function" on a doubled space. The basic ideas inherent in this picture 
were developed by Hormander (1971) and are also present in Miller (1974), while an 
elementary explanation of the geometrical situation is given by Littlejohn (1990). 

It turns out that this picture generalizes to linear maps between Hilbert spaces 
M : % — > H', where the two Hilbert spaces need not have the same dimensionality 
and where the map need not be unitary or even invcrtbible. The genralization involves 
again symplectic reduction, but in a different manner than that in which it appears 
in this paper. This situation arises in the comparison of two models of the 3j-symbol, 
which is a simpler version of the comparison between the 12 j- and 4j-models of the 




(93) 
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6j-symbol. The first is a 2j- or Clebsch-Gordan model, in which the Hilbert space is 
Cj 1 <£>Cj 2 , where the operator J 3 is a function of Ji and J2, that is, J 3 = Ji + J 2 . The 
second is a symmetrical or 3j-model, in which the Hilbert space is Cj t ® Cj 2 ® Cj 3 and 
the operator J3 is independent, but in which we are interested only in states satisfying 
the constraint Ji + J2 + J3 = 0. We will report on these investigations in the future. 

Roberts' (1999) derivation of the Ponzano-Regge formula effectively proceeded 
by writing the 6j-symbol as a scalar product (B\A), where \A) and \B) are given by 
Fig. [33J then writing the wave functions for \A) and \B) in the Bargmann (1962) or 
coherent state representation to obtain an integral representation for the 6j-symbol 
as an integral over C 24 . In the coherent state representation there is one copy of C 
for each degree of freedom, while in a 12j-model there are two degrees of freedom 
for each j, hence 24 degrees of freedom total. Roberts then used the stationary 
phase approximation to evaluate the integral. The Bargmann representation has been 
used in a similar manner in several recent asymptotic studies in the quantum gravity 
literature. 

In I and in this paper, however, we have mostly worked in a represenation- 
independent manner. Our emphasis on Lagrangian manifolds and other geometrical 
structures in phase space is part of this approach. From a certain point of view the 
coherent state representation is just another representation, albeit a complex one, so 
there is the question of whether it plays any privileged role or offers any advantages. 

For some purposes it certainly does, for example, the wave functions in the 
coherent state representation are polynomials that can be written down explicitly, 
and these in turn are useful for deriving generating functions and other things, as 
shown by Schwinger (1952) and Bargmann (1962). Another point is that in the case of 
SU(2), the coherent state representation appears naturally in the method of geometric 
quantization, that is, the holomorphic sections of Hermitian line bundles over the 
symplectic manifold S 2 are represented explicitly by Bargmann's entire analytic wave 
functions. 

On the other hand, it is not obvious to us that the identification of the 
stationary phase set or the other aspects of the semiclassical calculation are easier 
or more transparent in the coherent state representation than in the representation- 
independent approach of this paper. Moreover, if one focuses too narrowly on the 
stationary phase evaluation of an integral, one misses the geometrical structures that 
support the representation-independent approach. 

We had to assemble and extend several of the ideas presented in this paper for 
our derivation of the asymptotic form of the 9j-symbol, and we have studied other 
extensions such as "^-inserted" spin networks, in which a group element or D-matrix 
is inserted in the edges of the spin network. Such ^-inserted spin networks are basic 
amplitudes in loop quantum gravity (Rovelli 2004), taking one from the basis of spin 
connections to the spin network basis. The asymptotics of these amplitudes leads to 
piecewise flat manifolds similar to those introduced by Regge (1961). 

We will report on these and other developments in the future. 
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